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On  the  basis  of  the  modified  He  11 inger-Reissner  variational 
theorem  of  nonlinear  elasticity,  a  general  nonlinear  theory  is  devel- 
oped for  sandwich  plates  that  are  stressed  in  the  plastic  region. 
The  plasticity  equations  correspond  to  the  stress-strain  equations  of 
Hencky,  extended  to  large  deformations.   As  a  result,  the  material  is 
taken  to  be  isotropic,  and  the  loading  is  assumed  to  be  both  propor- 
tional and  active.  The  general  theory  is  then  extended  to  the  buckling 
of  sandwich  plates  by  applying  the  bifurcation  theory. 

In  addition,  by  omitting  certain  quadratic  terms  from  Green's 
strain  tensor,  a  simplified  large  deflection  theory  is  developed  for 
the  plastic  bending,  stretching  and  buckling  of  sandwich  plates. 
Furthermore,  the  theory  is  then  extended  to  those  sandwich  plates  which 
have  thin  outer  layers;  in  this  case,  the  Kirchhof f-Love  hypothesis 
is  applied.   Finally,  a  numerical  solution  is  obtained  for  the  plastic 
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buckling  of  a  simply  supported  sandwich  plate  with  very  thin  outer 
layers  and  a  soft  core  under  a  uniaxial  compressive  load. 
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1 .  .INTRODUCTION 

By  way  of  introduction,  it  probably  should  be  mentioned  first 
that  many  man-made  structures  consist  of  f lexural  members,  that  is, 
plates,  shells  and  thin  bars.   For  these  structures  to  behave  in  a 
predictable  manner,  the  ability  to  analyze  the  f lexural  members  is  of 
vital  importance.   In  general,  f lexural  members  fail  in  one  of  two 
ways,  either  they  become  overstressed  or  they  deform  excessively,  in 
particular,  the  buckling  phenomenon.  Often,  two  separate  analyses  are 
required  to  investigate  both  possibilities.  Well-established  branches 
of  mechanics  have  been  developed  which  deal  with  such  analyses,  for 
example,  strength  of  material,  theory  of  elasticity,  theory  of  plas- 
ticity, theory  of  plates,  theory  of  shells,  theory  of  elastic  stabil- 
ity, etc. 

In  this  work,  one  small  class  of  f lexural  members  will  be  con- 
sidered, namely,  a  three-layered  plate  which  is  commonly  called  a 
sandwich  plate.  Usually,  the  primary  purpose  of  a  sandwich  plate  is 
to  obtain  a  relatively  rigid  plate  that  weighs  less  than  a  solid  (one 
material)  plate  of  the  same  rigidity;  however,  sandwich  plates  have 
certainly  been  designed  for  other  purposes. 

The  sandwich  plate  has  found  numerous  applications  in  the  air- 
craft and  aerospace  industries;  as  a  result,  considerable  attention  and 
work  has  been  devoted  to  this  type  of  plate  in  the  last  twenty  years. 


Because  of  their  geometry,  it  is  natural  to  develop  a  theory  for  sandwich 
plates  by  extending  the  theory  of  ordinary  plates.   Unfortunately,  how- 
ever, there  are  many  theories  for  ordinary  plates;  for  example,  there 
are  theories  which  depend  upon  the  plate  geometry  (thick  or  thin),  and 
the  amount  of  deformation  (small  or  large  deflections),  and  the  state 
of  stress  (elastic  or  plastic).   Even  more  numerous  theories  can  be 
developed  for  sandwich  plates  because  of  its  structure;  for  example, 
consider  the  outer  two  layers  which  may  be  thick  or  thin,  or  different 
thicknesses,  or  different  materials,  etc. 

Because  of  this  multitude  of  theories,  there  comes  a  time  when 
it  is  desirable  to  develop  a  general  theory  which  encompasses  the  var- 
ious specialized  theories.  However,  by  their  nature,  the  equations  in 
a  general  theory  are  usually  too  complicated  to  be  solved,  and,  there- 
fore, recourse  must  be  taken  to  the  simpler  theories  for  solutions. 
One  of  the  advantages  of  reducing  a  general  theory  to  a  special  theory 
is  that  one  is  made  aware  of  the  factors  being  neglected.   In  contra- 
distinction, when  the  special  theory  is  developed  directly,  one  is 
usually  aware  of  only  those  factors  being  considered;  note  the  dif- 
ference.  Because  of  the  variety  of  methods  available,  a  general  the- 
ory can  be  developed  by  different  techniques,  for  example,  from  one  of 
the  variational  principles  or  directly  from  the  three-dimensional 
equations  of  nonlinear  elasticity.   Now,  the  fundamental  equations  of 
the  nonlinear  theory  of  elasticity  consist  of  the  following:   the 
equations  of  motion,  the  boundary  conditions  (both  stress  and  displace- 
ment) »  the  strain-displacement  equations,  and  the  stress-strain  equa- 
tions; these  equations  are  the  ingredients  of  a  general  plate  theory. 


Recently,  a  general  theory  for  sandwich  plates  was  given  by  Ebcioglu 
[l],  in  which  Hamilton's  principle  was  employed. 

Because  of  its  importance  to  this  work,  it  is  desirable  to 
discuss  some  of  the  salient  characteristics  of  Ebcioglu's  work.  The 
equations  were  derived  with  reference  to  the  undeformed  and  stress 
free  state  of  the  plate,  called  the  reference  state;  in  this  frame  of 
reference  (material  or  Lagrangian),  there  are  two  types  of  stress 
tensors,  the  Lagrangian  stress  tensor  and  the  Kirchhoff  stress  tensor. 
The  Kirchhoff  stress  tensor  has  the  advantage  of  being  a  symmetric 
tensor  and  was  used  by  Ebcioglu.   From  geometrical  considerations,  the 
most  general  strain-displacement  relations  are  the  Green's  strain  ten- 
sor (material  coordinates)  and  the  Almansi's  strain  tensor  (spatial 
or  Eulerian  coordinates),  both  of  which  are  nonlinear;  Ebcioglu  used 
Green's  strain  tensor.  From  a  physical  point  of  view,  Ebcioglu  took 
very  general  stress-strain  equations  which  include  thermal  strain  and 
are  applicable  to  elastic  anisotropic  materials;  the  equations  allow 
for  the  possibility  of  a  nonlinear  material. 

From  the  preceding  paragraph,  it  becomes  apparent  that  the 
general  theory  can  take  different  forms  deDending  upon  the  type  of 
stress  tensor  applied  and/or  the  frame  of  reference  used.  However,  in 
this  work,  an  alternate  form  of  the  general  theory  will  not  be  devel- 
oped; instead,  the  theory  will  be  developed  from  a  different  varia- 
tional principle,  namely,  the  modified  Hel linger-Reissner  theorem. 
The  modified  Hel 1 inger-Reissner  theorem  has  recently  been  employed  by 
Habip  [2]  to  develoD  a  general  theory  for  ordinary  plates,  and,  later, 
Huang  [3]  applied  the  theory  to  a  sandwich  shell. 


It  should  be  mentioned   that   for   a  general   theory  there   are 
no  uniqueness   theorems  which  state     that  the   solution    is   the  only  pos- 
sible  solution;    this    is   as    it  should  be  because    it    is   well   known  that 
there   are  different  deformations   possible   for  a  given  environment , 
commonly  known  as   the  buckling   phenomenon.     Therefore,    a  special   for- 
mulation  is   required   to  predict   the   onset   of   buckling.     To   the  •author's 
knowledge,   the  general   theory  for  sandwich   plates   has   not   yet  been 
extended   to  the  buckling   problem;    the  buckling  equations  will   be   devel- 
oped   in  this  work. 

Although  the   stress-strain  equations   employed  by  Ebcioglu  may 
be  nonlinear,    it    is  well   known   that   for  most   structural  materials    (i.e., 
steel,    aluminum,   etc.)   a  nonlinear  stress-strain  equation   implies   that 
the  material  has   been   stressed    into  the   plastic  range.      Because   the 
theories   of   plasticity  and  elasticity  are   not   the   same,   different  gen- 
eral  theories   should  be  developed  which  are   applicable  to  both   types   of 
plate   problems.      In  this  work,   we   are   interested   in  developing   a  gen- 
eral  theory  which   is   based   on  the  equations    of   plasticity. 

In  the   field   of   plasticity,   there   are  many  theories,   usually 
classified   as   the   incremental    (or   flow)   theories   and  the   deformation 
(or   total   strain)   theories.     The  book  by  Hill    [U ]  deals   almost  exclu- 
sively with  the   incremental   theories    (in  particular,   the  Reuss   equa- 
tions)  because   these  theories   describe   the   actual   behavior   of  materials 
better   than  the   deformation  theories.     However,    the  deformation   theories 
are   simpler  mathematically  and,    therefore,    are   ofted  used,   particularly 
Hencky's   equations.      It    is  well   known  that  Hencky's   equations   can  be 
derived   as   a  special   case   of   Reuss*    equations,   namely,   when  the  principal 


axes   of   stress   do  not  change   and  the  stress   ratios   remain  constant, 
see   Smith   and   Sidebottom    [5];    this    is    a   practical    condition   that  occurs 
quite   often   in  plates.      In  Hill's    [4]   comments   on  Hencky's   equations, 
he  mentions    that   the   equations   have   not   been  extended    to  the   case   of 
large   strains    (e.g.,   Green's   strain  tensor);   however-,   as  we   shall   see, 
the   foundations   have  been  established   for   large   strain   theory. 

The  modified  Hellinger-Reissner  theorem  assumes   the  existence 
of   a  strain  energy  density  from  which  an  explicit  form   for   the  stress- 
strain  equations    is    obtained.      However,    a   strain  energy  density  implies 
that   the  material    is   elastic,   not   plastic;   but,    according  to  Novozhilov 
[6],    Kachanov  has    proved   that   for   an  active  deformation   an  elastic- 
plastic  body   is    indistinguishable   from  an    ideally  elastic   body,    both 
with  the   same   stress-strain  curve.     Therefore,   the  concent   of   a  strain- 
energy  density  can  be  employed   to  a   restrictive   class   of   plasticity 
problems.      Furthermore,    for    large   strains    (Green's   strain   tensor), 
Novozhilov    [6]   has    shown   that   equations    corresponding   to  Hencky's   equa- 
tions  can  be   obtained   in   terms   of   a  strain  energy  density.     The   proce- 
dure  of   Novozhilov  will   be  used   in  this  work;    and,    furthermore,   explicit 
expressions   for   the   strain  energy  density  and,   hence,   Hencky's   equa- 
tions will   be  derived. 

The   theory  for   the   plastic  buckling  of   plates    is    largely  a 
result   of   experience   gained   in   the  study  of   the   plastic  buckling  of 
columns.      As    is   well   known,    two  different   philosophies   had   prevailed 
for  many  years   concerning   the   plastic   buckling   of   columns,   which  were 
based   on   the   assumption  of  whether   or  not   strain  reversal    occurs   during 
buckling.     The  two  different   approaches  were  unified  by  Shanley   [7] 


in  1947.  As  discussed  by  Gerard  [8],  von  Karman  concluded  from 
Shanley's  work  that  the  stability  limit  should  be  redefined  for  plastic 
buckling  as  "the  smallest  value  of  the  axial  load  at  which  bifurcation 
of  the  equilibrium  positions  can  occur,  regardless  of  whether  or  not 
the  transition  to  the  bent  position  requires  an  increase  of  the  axial 
load."  Therefore,  the  bifurcation  theory  for  buckling  will  be  employed 
in  this  work;  a  formal  statement  of  the  theory  is  given  in  Sect.  3.4. 
It  is  important  to  note  that  in  the  bifurcation  theory  the  loading  is 
active. 

Based  on  the  above  discussion,  a  general  theory  for  sandwich 
plates  which  includes  buckling  is  developed.  The  theory  is  based  on 
stress-strain  equations  of  the  Hencky  type,  extended  for  large  strains. 
The  theory  is  developed  in  Sects.  2  and  3  of  this  work.   In  addition, 
some  specialized  formulations  are  considered. 

In  bending  problems,  the  lateral  displacements  are  usually 

much  larger  than  the  displacement  components  in  the  plane  of  the  plate. 

* 
On  physical  grounds,  it  has  been  argued  that  often  certain  quadratic 

terms  associated  with  the  inplane  displacements  may  be  omitted  from  the 

strain  tensor.   In  Sect.  4,  a  modified  theory  is  developed  in  which 

only  those  quadratic  terms  of  the  lateral  displacement  are  retained  in 

the  strain  tensor. 

Because  most  sandwich  plates  are  made  with  thin  outer  layers, 

special  consideration  is  given  to  this  type  of  structure.  As  mentioned 

by  Fung  [9],  probably  the  most  important  discovery  in  plate  theory  was 


*  r  i 

See  Novozhtlov  [6  J. 


Kirchhoff's   first  assumption   (stated   in  Sect.    5)  which  reduced  the 
equations    to  a  form  where   practical   solutions   could  readily  be   obtained. 
Kirchhoff's   assumption   along  with   Love's   assumption  are  basic  hypotheses 
in  current   thin   plate   theory.      In  Sect.    5,    a  theory  for  sandwich  plates 
with  thin  facings    is   developed  which  utilizes   the  Kirchhof f-Love  hypo- 
thesis   along  with  the  modified  strain   tensor  described   in  the   preceding 
paragraph. 

Finally,    in  Sect.    6,    a  sample   problem   is  worked   out    in  detail 
and   compared  with  the  work   of   previous    investigators. 


2.      THE   VARIATIONAL   PRINCIPLE 

It   is  well   known  that    in  the   nonlinear  theory  of  elasticity 

the  basic  equations   can  be   obtained   from  variational    principles,    in 

* 
particular,   Reissner's   theorem.      Because   of   the  early  work  of  Hellinger 

[10],    the  variational   principle   is   often   referred   to  as   the  Hellinger- 

Reissner   theorem,    see  Truesdell   and  Toupin    [12 ].        Rather   than  work 

with  the   theorem  as   stated   by  Reissner    [13],    an  alternate   form  developed 

by  Washizu    [lU]    is  more   convenient   for   our   purnoses,    referred  to  as   the 

modified   Hel 1 inger-Reissner   theorem.      Habip    [2]   has   used    the  modified 

Hellinger-Reissner   theorem  to  develop  a  general   nonlinear   plate   theory; 

in  this  work,    the   procedure  used  will   be    quite   similar   to  that   given 

by  Habip.     The  variational   principle   as   given  by  Habip    is 

"the  modified  Hellinger-Reissner  theorem   asserts   that   the  varia- 

tional  principle 


(\{f-fl)lvi  M 


=   h 


(s'J'rt.-z')<tr 


{S'j(utl.  +   Vjli  +    Vrl:   Vr\:)<LV 


(2.1) 


*  r         i 

According   to  Naghdi   [11],   Hellinger's  work  did  not    include 
the  boundary   integrals. 
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where  y.  .,  S   ,  v  and  S+  are  varied  independently,  is  equiv- 
alent to  Cauchy's  first  law  in  V,  to  the  stress  boundary  condi- 
tion on  the  part  A  of  the  boundary,  to  the  displacement  boundary 

condition  on  the  remaining  part  A  ,  and  to  the  stress-strain  and 

o  v 

strain-displacement  relations    in     V,   when  the   symmetries   of  y.  . 

and  S   are  both  used." 
All  of  the  various  quantities  in  the  above  theorem  are  referenced  to 
a  body  that  is  unde formed  and  stress  free,  known  as  the  reference 
state;  for  example,  the  symbol  V  denotes  the  volume  of  the  undeformed 
body,  i.e.,  the  volume  in  the  reference  state.  With  the  understanding, 
hereafter,  that  all  of  the  quantities  are  with  respect  to  the  refer- 
ence state,  the  meaning  of  the  various  symbols  used  in  the  statement 
of  the  theorem  is  as  follows: 

A   is  the  portion  of  the  boundary  of  the  body  upon  which  the 
displacement  vectors  are  prescribed  (given); 

A   is  the  remaining  part  of  the  boundary  upon  which  the 
stress  vectors  are  prescribed; 


o 


f   are  components  of  acceleration; 

F   are  components  of  body  force  per  unit  mass; 


o 

S    are  contravariant  components  of  the  stress  tensor; 

S+   are  contravariant  components  of  the  stress  vector; 

v    (v*)  are  contravariant  (covariant)  components  of  the 
displacement  vector; 

V   is  the  volume  of  the  body; 


are  covariant  components  of  the  strain  tensor; 
6    is  the  variational  symbol  used  in  calculus  of  variations; 


Vij 


0  is   the  mass   density  of   the  body; 

£  is   the   strain  energy  density  per  unit  volume; 
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is  the  covarLant  derivative  symbol  used  in  tensor  calculus; 

is  the  tilde  symbol  which  is  used  to  denote  prescribed 
quantities . 

When  the  indicated  variation  is  performed,  Eq.  (2.1)  can  be  put 
into  the  form 


{slJ(S^vy}\.^J0Fr-r)}6vr  JLV 
4     (v.-  v.  )   hs[  dA 


where  Gauss'  theorem 


y 


(2.2) 


(SiJ  6v;)l.  dV  =      SiJ 'Sir-  m.  J.  A  (2.3) 


has  been  employed  and  where  the  strain  energy  density  has  been  assumed 
to  be  a  function  of  the  strain  tensor,  i.e., 


£   =  Z(*ti). 


J 

In  Eq.  (2.2),  6^  is  the  Kronecker  delta,  and  n.  is  the  external  unit 
normal  vector  of  the  boundary  surface.  Then  from  the  generalization 
and  extension  of  the  basic  lemma  of  calculus  of  variations  (hereafter, 
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referred  to  as  the  fundamental  theorem  of  calculus  of  variations),  as 
indicated  by  Weinstock  [15],  Eq.  (2.2)  yields  the  fundamental  equations 
of  the  nonlinear  theory  of  elasticity  in  terms  of  the  reference  state 
as  follows: 


the  equations  of  motion  in  V  (Cauchy's  first  law), 


J*/ ,  i 


#{?k*v\)\  yAr-D-o-, 


the  boundary  conditions    on     A    , 

os 


<>>/!  i 


Stt-S>iS'*(Vk*V\)      -~    0   ;  (2.5) 


the  boundary  conditions   on     A    , 

1  O    V 


V-    -     V-    =    0;  (2.6) 

the   strain-displacement   relations    in     V, 

o 

Y..  =J{VJJ   -  Vjir+Vl;   Vrlj)  ;  C2.7, 

the  stress-strain  relations    (in  terms   of   the   strain  energy  density) 
in     V, 


6  2      hi-  )>Y:.)        ' 


2l^i      *V 


* 

These   equations    are   geometrically  nonlinear;    and   depending   uoon  T   , 
the   equations   can  also  be   physically  nonlinear.      Later,    in   Sect.    3.3.4, 
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it  will  be  shown  that  by  a  proper  choice  of  T,     the  equations  are  appli- 
cable to  a  certain  class  of  plasticity  problems.   Perhaps  it  should  be 


or 


mentioned   that    in  the   literature     y.  .      is   known  as   Green's   strain  tens 

and     S.  .      is   known   as   Kirchhoff's   stress   tensor. 
i-j 

Examination  of   the  modified  Hel Iinger-Rei3sner   theorem  shows 
that  by  certain  modifications,   which  are   practically  self-evident,   equa- 
tions  corresponding  to  alternate   theories   can  be   obtained.      For  example, 
in  the   classical    theory   of   elasticity,   the   strain  displacement  equation 
is 


Xj =  iivtij+vjiji 


it  is  obvious  that  this  expression  can  be  obtained  from  the  theorem  by 
omitting  the  quadratic  term 


vrl;  ur\ 


from  the  third  integral  in  Eq.  (2.1).  Thus,  Eq.  (2.1)  would  read 


pAf-fl^iv-  h\\{s%-t)dv 


Ic'J 


Is    ^lj  +  Vjh)dV  +\a  S>,  dA 
+L    Si  (  V,  -  vt  )  dA 

which  after  carrying  out  the  indicated  variation  becomes 
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( K  -  S lif,j  )SVidA+\    (vr  *i  )«i  dA 


%  -{(Vi/.*  vjUnSs'JJV 


The  application  of  the  fundamental  theorem  of  calculus  of  variation, 
indeed,  yields  the  classical  equations  of  elasticity,  provided  of 
course,  that  S    is  now  interpreted  as  the  Eulerian  stress  tensor 
and  T.     is  chosen  to  yield  Hooke's  law. 

In  the  study  of  plates,  theories  have  been  developed  which  are 
based  on  both  the  general  nonlinear  equations  and  the  classical  equa- 
tions.  In  addition,  intermediate  theories  have  been  developed,  for 
example,  the  work  of  Ebcioglu  [16]  in  which  a  strain  tensor  was 
employed  as  follows: 

Kj  =  I  (Vi/j  +  vj/i  -  n,i  irs,j) 

where   the   comma  denotes    the   partial    derivative;    a  more   restrictive 
form  of    the   above   strain   tensor    is   used    in   the  well-known   von  Karman 
theory.      Because  virtually  all    plate  theories   are   based   on   simplify- 
ing assumptions,    it   is   sometimes   difficult   to  establish   a   consistent 
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set  of  equations,  particularly  when  the  equations  are  developed  by  the 

* 

direct  method     which   is   often  employed.      In  general,    this   difficulty   is 

overcome   by  using   the  modified  He  11 inger-Reissner   theorem  because  by 
introducing   the   simplifying  assumptions    in  a   systematic  manner   the 
theorem  automatically   produces    a   consistent   set    of   equations. 


In  the  direct  method   the  basic  equations   are   obtained  by 
considering   the   stress   resultants   acting  on   a   plate  element   and   then 
applying   the  equations   of   equilibrium. 


3.   GENERAL  NONLINEAR  THEORY 
OF  A  SANDWICH  PLATE 


3.1.   Description  of  the  Sandwich  Plate 

In  this  work,  a  sandwich  plate  is  a  three-layer,  laminar 
structure  of  constant  thickness.  The  layers  may  be  dissimilar  mater- 
ials of  different  thicknesses;  however,  the  thickness  of-each  individ- 
ual layer  is  constant.   In  addition,  each  layer  is  intimately  fixed  in 
relation  to  the  other  so  that  the  sandwich  plate,  as  a  whole,  is  a  con- 
tinuous medium;  thus,  the  displacements  are  continuous  throughout, 

including  across  the  interfaces.  As  a  result  of  the  application  of 

* 
Gauss'  theorem  to  obtain  Eq.  (2.2),  the  plate  must  be  a  regular  region; 

therefore,  there  are  no  cavities  or  holes  through  the  plate.  Also,  it 
is  desirable  to  require  that  the  edge  surface  is  a  portion  of  a  right 
cylindrical  surface.   Naturally,  the  terminology  "plate"  implies  a 
flat  structure  (as  distinguished  from  a  shell  which  has  a  curved  sur- 
face^  the  thickness  of  which  is  small  as  compared  with  its  other  dimen- 
sions. 

** 
To  properly  describe  the  sandwich  plate,  a  coordinate  system 

must  be  chosen;  naturally,  curvilinear  coordinates  should  be  chosen  to 

utilize  the  general  form  of  the  variational  principle.  However,  to 


*  r    -i 

See  Kellogg  [17  J. 


** 

All   coordinate  systems  will   be   taken   as   right-handed, 


15 


16 


take  advantage  o£  the  geometry  of  a  plate,  the  curvilinear  coordinates 
can  be  specialized  without  introducing  a  loss  in  generality  to  the 
theory;  these  specialized  coordinates  and  the  rectangular  Cartesian 
coordinates  to  which  they  are  referred  are  shown  in  Fig.  1. 


Fig.  1.   Plate  coordinates 


Here,  the  9  coordinate  curve  is  rectilinear  and  orthogonal  to  the 

1       2 
8  and   9  coordinate  curves;  furthermore,  the  functional  relationships 

between  the   9  and  x  coordinates  are 


e   =  d"(xP) 


n3   -    3 

U       —  X      -h   const. , 


and 


/*=    S*(*P< 


X      (  6     )   -      X      -     6      +   const., 
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where  the  convention  is  now  introduced  that  the  Greek  indices  (iv,  P,v, . . .  ) 
take  only  the  values  1,2.   Note  that  all  planes  parallel  to  the  upper 
and  lower  faces  of  the  plate  are  described  by  equations  of  the  type 

U        "  const. 

By   introducing   a  translation  of   coordinates,    any  plane   parallel   to  the 
faces    can  be  written   as 


-2  3 

S       -     8     -     const.   =   0.  (3.1) 


In  plate  theory,  the  usual  procedure  is  to  write  the  equations  in  terms 
of  coordinates  such  that  the  middle  plane  is  described  by  an  equation 
like  Eq.  (3.1).   Therefore,  it  is  convenient  to  choose  the  coordinate 
system  for  the  sandwich  plate  as  shown  in  Fig.  2  where  the  middle  planes 
of  the  upper,  middle,  and  lower  layers  are,  respectively,  described  by 
the  equations 


»*3  -3  „3 

6=0,  0=0,  6=0;  (3.2) 

and  are  related  by  the  equations 


'e3  =  e  -  £(h  +  'h), 


i*3  n3         i 

(3.3) 
"~3  ^3 


e=d+±(h  +  %) 


where   'h,  h,  "h  are  the  thicknesses  of  the  upper,  middle,  and  lower 
layers,  respectively.  Hereafter,  the  notation  of  a  prime,  a  bar,  and 


X 


18 


-P~X 


Fig.  2.   Coordinates  for  sandwich  plate 
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a  double  prime  will  be  used  to  distinguish  quantities  associated  with 
the  upper,  middle,  and  lower  layers,  respectively. 

For  the  above  coordinate  system,  the  metric  form 


4>  =  cjr  do  do 


becomes 


,,/3      .     /   /„3t* 


§  =  9w/9  je*je*  +  ue°y 


where,  from  the  definition  of  the  metric  tensor, 


-  L    \p>« 


Jxr   Jx 


%P~  fa  her  he"  9 


a 


*3 


=    0 


% 


33 


1.0  . 


(3.4) 


As  a  result  it  is  easy  to  show  that  all  the  Christoffel  symbols  with 
the  index  3  become  zero,  i.e., 


id, 


i 


J 


r  3  1 


0 


(3.5) 


J3r\3Ji =  vj(  =  ° 


also 


0 


*&  C&  °<3  33 


(3.6) 
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Because  of  the  above  characteristics, 


<*  °( 


VLl3--      V\3      ,  Vil3     =     V;,3     ,        (3.7) 

v'h  =  v3,i    »        v3li  '  v3ii  , 

where  the  comma  denotes  the  ordinary  partial  derivative  with  respect 
to  the  9  coordinates. 

The  coordinate  system  is  the  Euclidean  space  case  of  the 
geodesic  normal  coordinate  system  discussed  by  Synge  and  Schild  [18 ]; 

as  a  result,  when  a  space  tensor,  say  S   ,  is  split  up  into  the 

^oB  „o3   „33   .  * 

components  S   ,  S   ,  S   ,  the  components  are  known  as  subtensors , 

subvectors  and  subinvar iants ,  respectively. 

3.2.   Displacement  Functions 

Throughout  the  region  of  the  plate,  the  displacement  is  assumed 
to  be  a  sufficiently  well-behaved  function  that  can  be  expanded  in  a 
Taylor  series,  at  least  to  the  first  two  terms.   For  our  purposes,  it 
is  desirable  to  expand  v  in  terms  of  9  about  the  plane   9  =  0,  i.e., 

v'(dJ)=   vL(d",o)  +d3[vi(6\o)     ]. 


* 
Some  authors  use  the  terminology  surface  tensors,  e.g., 

see  Green  and  Zerna  [19]. 


M 
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By  introducing  the  notation 

<il  iff*)  =    vL  (&",  d) 

the   displacement   can  be  vnritten  as 

Vl(eJ)  -    nl '  (?)  +  83fl(6*).  <3..) 

The  function  u  represents  the  displacement  of  a  particle  in  the 
middle  surface.  As  mentioned  by  Novozhilov  [6],  the  function   * 
characterizes  the  direction  of  a  fiber  in  the  strained  state  which 
was  initially  normal  to  the  middle  plane.   Note  that  fibers  which  were 
initially  normal  to  the  middle  plane  are  not  necessarily  normal  to  the 
middle  surface  after  deformation,  as  assumed  in  many  plate  theories. 
By  extending  the  above  to  a  sandwich  plate,  the  displacement 
of  particles  in  each  of  the  three  layers  can  be  written  as 


'vc('eJ)  =  'uc ($"■)  +  'e3  f^e")  , 

vc  (§J)  =  u'(e")  *  e3  r'(#V,         o.« 
"vL("eJ)  ='tti'(f)+  Vy  (6*). 

By  use  of  Eqs.  (3.3),  the  equations  can  be  rewritten  in  terms  of  the 
coordinates  of  the  middle  layer,  i.e., 
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'v£(3J)=  ki(**)+[9-{(*  +  'l')]'ri'(f) 


~3     -L,«« 


VL  (6J)  =    LLL  (0   )    *    0      ?'(*"),  O.io) 


_  3 


v<  ( 9j)  =  v  c<r> [$+i(**  '*>)]  Y  (£>*)  ■ 

These  equations  can  be  related  further  by  making  use  of  the  fact  that 
the  displacements  must  be  equal  at  the  interfaces  (continuity  condi- 
tions) 


'vl 

= 

*?<', 

at 

V 

= 

*■', 

at 

-»'  -  - 1 

(3.11) 


Substitution  of  Eqs.  (3.10)  into  (3.11)  gives 


V=  £ 


+  2  I  h  if    +  h  w    J 


(3.12) 


v  =  v  -  ±[h  yL  +  "h  y] 


or 


ipl=U'u}-uL-{('h'y  +  "hy)] 

Therefore,  the  displacements  can  be  written  as  either 


(3.13) 
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VL=     LLL+    \    (lpl-    '(f/L  )  +  83  if1 


I 
TfC  =      HL    +     6~    ifj'  (3.14) 


Z-A  rrc     .      E3     -l 


or 


r  ? 

If      =     LL     -    ~f    (  (jJL  -    tfJ    )    +    6      lfJL 


'vL  =  u1'  +  [33-{(h  +  'h)]  'fL 

-3  (3.15) 

h 


In  applying  the  variational  principle,  either  Eqs.  (3.14)  or 
Eqs.  (3.15)  can  be  substituted  into  Eq.  (2.2).   If  the  former  are 
used,  the  resulting  exDressions  will  be  in  terms  of  the  quantities 


-l  -    t  It  a.      L 

or,  if  Eqs.  (3.15)  are  employed,  the  results  will  be  in  terms  of 


I       I  il       l  I        i  a       L 

ll  ,   a  7   if1  ^f  . 

In  either  case,  there  will  be  twelve  unknown  displacement  functions 
to  be  determined.   Regardless  of  which  set  of  variables  one  chooses 


2U 


to  work  with,  the  resulting  equations  can  be  reformulated  in  terms  of 
the  other  variables  through  Eqs.  (3.12)  or  (3.13);  this  will  be  demon- 
strated later  for  a  simple  case.   Because  Eqs.  (3.1U)  are  more  compact 
than  Eq.  (3.15),  the  equations  will  be  formulated  in  terms  of  the 
variables 


—  L  —    L  I         I  ll         C 

*     ,     V     y      f     7     f 


3.3.   Fundamental  Equations 

3.3.1.   Equations  of  motion 

Because  the  variational  principle  leads  to  a  rather  lengthy 
equation,  i.e.,  Eq.  (2.2),  which  is  awkward  to  handle  in  its  entirety, 
it  is  desirable  to  split  ud  the  equation  if  possible.   Examination  of 
Eqs.  (2.2)  and  (2  .U)   shows  that  the  equation  of  motion  is  a  direct 
consequence  of  the  first  integral  of  Eq.  (2.2).  Therefore,  the  equa- 
tion of  motion  can  be  obtained  from  the  equation 

y  J   '  (3.i6 

When  written  in  terms  of  the  subtensors,  as  discussed  in  Sect.  3.1, 
Eq.  (3.16)  becomes 
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r 


W, 


.3      n3 


(3.17) 


For  future  use,  it  is  convenient  to  rewrite  Eq.  (3.17)  as  follows 


(3.18) 


where 


1/8  _   (St 


P3, 


*p-^(^V'lt)  +  s<"(^v\3) 


3i3-s'3(&i  +  v%)+S33(&i.V\3). 


(3.19) 


The  total  volume   V  of  the  sandwich  plate  consists  of  the  sum  of  the 
o 

upper  layer  volume  'v,  the  middle  layer  volume   V,  and  the  lower  layer 

o  o 

volume  "V.   Therefore,  for  a  sandwich  plate,  Eq.  (3.18)  becomes 
o 
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'<y  (3.20) 


Note   that 


d'\/  =  d'e3dA,  civ-  d$3J/iy  ft/-  d"d3oiA    (3.2D 


where 


Then,  with  Eqs.  (3.21),  the  integrals  of  Eq.  (3.20)  can  be  combined 
to  obtain  the  equation 


.a 


■c 


(3.22) 


where    the    limits    are 


a-  -  *  +  2    '  A       2   ' 

(3.23) 

c  =  -  I  ,  <•  -  -  "h  -  I  . 
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Expressions  for  the  displacements  were  established  in  Sect.  3.2;  thus, 
substitution  of  Eqs.  (3.14)  into  Eqs .  (3.22)  gives 


a 


h  ' 


LL 


vt 


■c 


+ 


Id 


(e'+i)ib«L  +  '"l  if 


(3.24) 


b 


Li 


3 


*rim+-]+J  *'[s*l+~  i 


L'n---pr3+C(63-kn--pk 


A  f"r\3'1 
2 


-3 
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where 


W-s"Lii  +  u%  +  i(r%-y/f) 


_-5  / 


+  g  if,*/,]  t!sfi3(&i  +  'f< 


o'P.  §<>* 


6<   -  u%  +  §  ft  ) 


-  133 


+  &3yit]  +  V3(si  +  y 


k"  =  S«[6i  +  u%  +  !(?%-?%.) 


■=.  13 


E 


+  e3'p%]+  's33(bi  +  ?■), 


-  33 


(3.25) 
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Note   that 


■  a 


jb 


'£*\3  Je> 


'£<3k 


b    > 


a 


;£ 


£°%  03JB3* 


(*'k«t  - 


-a  (3.26) 

/    <*3     ,  -3 

E      dd   . 


Therefore,    Eq.    (3.24)   becomes 


,-J 


f> 


c    -  i-  ->d '  p 

h 


+  ik*')-b+(£**Je  + 


(*"%}* 


UL 


■b 


-h 


Jb 
"C 


h[b^--:]js3+r§3[D\--qjs3 


y 


I 1  o\ 


b 


l0o«L+~]j83-[£*s, 


-*3    ,  -3       £    ,i  jj 


6    +  -IE 


I 


a 


+ 


+ 


(9aE"£~U*'"£}i 


?* 


"( §-§)[  tf% +-]dd3-f'£«3S3 


3  *o(3  \a  h    /  'c/3 


a 


+  (i'Z"\  -iC^Ttl  h'n 


-3 


hj  (8-*±)[b%*-]JS-j  r3M   + 
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_5 


d3 


+    (6     £'  'I 


a   *$Ce"£}ih 


'3  in3 


'b' 

■b 


'fAr-mjs 


+ 


+ 


F33  f     +     f    F33^  /  1-33  \C     \      r   - 


h  r  u3<* 


-6 


+i  Ubls-1" rW  s'[d7s..-]Jg 


Jhn-]^3  -l^33^3 


z 


b 


c 


j.  t  /  'rr33f  +  /  a3  F33  \L  /  "    33 


hf 


3 


a 


+ 


(9 


7\3_    h    \[  'n3Mi 

'    2 


a 


,b.      <.J~ls~]JS'-J    'e"j§3 


+(i'*°):-§('£-):}m 


+ 


33, u  h     ,i     33  (X 

3     -  rC 

-s     l    .  r u 


>J 


e  +  IW/s-JJ/'-lk'sg3 


He3'E33p{(E%}yf3]d/i-o 


(3.27) 
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Next,  stress  resultants  are  defined  as  follows 


a  fb  rC 

.gp       §<fi    hyp 


I 

-3 


-3 


6       (JO      ■      (3.28) 


rJ,2 


mi 


V,  t\YV 


Jb        -c 


d     J 


-5 


(3.29) 


Jh  Jc  J  J  ) 


(3.30) 


r  '  •  -\ 

a.    3   a   ,   cl 


.a 


-A 


^1/^UHri^ 


r       c       _  £     //      / 

or  ?  us  ^  us- 


1       I  I       I  O    ~  '     I  I  "    »  J    I 

[Jh  JC  '  a 

-a  rb  rc 


_3 
(3.31) 


r 


bfy\fj'i]d 


i 


b  F\\    -pFW    "p  rL)\_i)JOi 


& 


(3.32) 


"£i3,S3'£iS)     ■      Vdy 
p\mCi\={EL\93£L3\,  *t\8%(-iJ1  l\-    (3.33, 


(        I      L  I      t    \ 


m 


Z 


\    j 


"k 

^  2 


•   L   .      !*<■        l^i  ',  ,rL  '  rL      .    'nL    \ 


P 


c 


a  ^  i      »lni       " , ,  r-  £*  "r> L       "n L 


)    .         (3.34) 


1,  t 


b  + 


y-,y,  -ur+zc-;c) 
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By  applying  the  fundamental  theorem  of  calculus  of  variations  to 

Eq.  (3.27),  and  by  using  Eqs.  (3.28)  -  (3.3*0,  the  equations  of  motion 

are  finally  obtained,  and  can  be  put  in  the  following  form: 

{(Wr+  N?r+W)(8*  ♦2*7, )  ♦  (/^r-  \  tf*)Ylt 


Id  —  o(  ii    o( 

a,   -r-  a    +  a. 


'■)  (3.35) 


{[  I  ('«"-  W<)  ♦  *<»]($  ♦  u\)  ♦  *>*  |  ^)H 


./A 


v 


2  ^         2       /r  if 

J  {CL    ~  &    J   +   /W        •  (3.36) 
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•MP*-EW)(g<t+B. 


°<l     +    h     -Aj      \   _    h    ^/3  Lrt 


ir  +  'ifit)'      9r'f 


■+ 


l(ifW-hW+'**fy%}L  +  CT'ir'„uyf 


'*  r  <* 


-<?[t>z1f+Ur%-y/*)]-ir'+'c 


2  cl  -t  /m      } 


(3.37) 


(w+iwxiun-$F/rj*s<n> 


41    U  h    ",J  ».M 


h  \2 « 


+[(i;  w, iwu -k^yuW  <t% - -n»)y 


h   *  «!       "  °< 


-o'[s;  +  u%-l(f%-'f%)]  +  ly*'c 


(3.38) 


(>#*+M"+W)a3t,+  (W-iWr,,, 


+ 


P3  +p3   +  p3     =      ^3^   ***   +   &, 


(3.39) 
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'[(iftfoW)**"]?,,  -  lew*  iwyf3tr 


t 


h  f/s 


z 


tfo+h)- 1  <?(»"&% +  (T%.-n")o+h) 


(3.40) 


XnKr)lf,*Crf'<fX<+k) 


h    i  > 

2    Cc^,    +   /W 


3      ? 


(3.41) 


AfP+k'A/fi')/,7  -h 


+ 


in    "  <> 

-zCL3-t    Aft3     . 


(3.42) 
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3.3.2.   Boundary  conditions  , 

It  is  obvious  that  the  boundary  conditions,  i.e.,  Bqs.  (2.5) 
and  (2.6),  are  a  direct  result  of  the  second  and  third  integrals, 
resoectively,  of  Eq.  (2.2).  Therefore,  as  in  the  preceding  section, 
the  boundary  conditions  can  be  obtained  from  the  equations 


(3.43) 


(3.44) 

The  total  boundary  A  of  the  sandwich  plate  consists  of  the  top  sur- 
face ^A,  bottom  surface  "A,  and  the  edge  surface  ^.   In  addition,  the 
edge  surface  consists  of  portions  of  the  upper,  middle,  and  lower 
layers;  thus, 


/    =7/f  V  +-(/+e"*/'1)-  <3-"5> 


Let  us  first  establish  the  stress  boundary  conditions  which  are 
obtained  from  Eq.  (3.43).   In  terms  of  subtensors,  Eq.  (3.43)  can  be 
written  as 
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or 


-3//   ,*  _«•/     ,  33       j        I    , 

~S>3  SUV3*  -J»3  $*0*  V3,3)]^}  'JA-O, 

J  A  ' 

0     5  (3.46) 

+  (Si  -j,    /**  -^3B33)bV3]dA  =  O 


where  Eqs.  (3.19)  have  been  employed.  Because  of  the  geometry  and 
coordinate  system,  observe  that: 


-  A . 

m„     -    0  .  /Y),    -     I.  O   : 

-  -A  > 

fa  :  °  ,      fa  '  '  I-  ° ; 


and  on  _/n  , 


j 


/3    *  /3     c      p     1         o       3        o     3        °      3 
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Also,  note  that  the  edge  integrals  can  be  written  as 


a 


)  clA    =J  J    {■■■)  J6SJ^  i 


■h 


A,  .-U  Jd 


£    5 


Therefore,  Eq.  (3.46)  can  now  be  written  as 


0    6 

(3.47) 

Next,  by  substituting  the  displacement  functions  as  given  by  Eqs.  (3.14), 
Eq.  (3.47)  becomes 
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'A 


'3i-'EM')[hs.siir4  +  (ss-k)&h] 


*  (Si  -  'E")[&U3  +  §$f3+(6>-l)S  fa]}  JA 


'Q. 


n3       h\/'?4 


+Jh(*-l)(s:-<*/o*)6'r< 


c 

-5 


-76  ^c  Jd 


tL 


jb  yJ~  z)(S*  ~^>  'D     )  tfa 

tf  (A  f )( J,3  -JS  V?)  I'yAs'j*    =  0 

d  J  (3.48) 

where  Eqs.  (3.25)  have  been  employed. 
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Now,    let   us    introduce   the   following   definitions: 


'   >'L  >n^ 


"p l     "c  l 


S*      5       &     S* 


I  \ 


-I 


3  i>       '    i 
-      ~  C    (   0 


6    *      ■)        ®      S  4 


a 

AT  6   =  <  J 

b 

V  ) 


(3.49) 


{  r  r  r 


.a 


>=< 


^  (     \ 

rb  ■    rc        / 

-6         Jc  ^ 


(9 


Then,   with   the   above   equations    and   Eqs .    (3.28)    -    (3.34),    the   boundary 
conditions    are  obtained    from  Eq.    (3.48)    by  applying    the   fundamental 
theorem  of   calculus    of   variations.      The   stress   boundary  conditions 
are  as    follows: 

on   the   top  surface, 


I      L  l     I 


-  0 


I     L  It 

2C*     -   f       =    °i 


(3.50) 


on  the   bottom  surface, 


,P*     +  ,P        =     0  7 


II      I-         ,        II     L 

A  +  ,c   -  o 


(3.51) 


on   the  edge   surface, 
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'  -  -        "X  *       a    < 


$  f'*tP.  r,P*  \,fi< 


t   L 


s   +o$  +  05     =  m   \  (  N'  *N  v  *  A/1'"  )  (  S  t  +  H * 


°s 


+ 1 4  {^^hM^f^  +  (m**-  {W*)Ylx 


Z 


r 


-t 


t*+  b/i/fi*)  Ml    -h'^  in*      7^ 


(3.52) 


4  /  ^-   C*  )  +    /*  =   ~    J  (A  I  *A/ fi*    U*  ^A7^H/<<        /7»S 


l(j-r)^=MLI(^-^)+^](K+ 


LL 


/ 


+ 1 '  *j  (y'f " 


2 


n^L  +  kC/^-h'^v,,^ 


ttrjflt+KW-r/ryflt 


mP+  f  ^9  yi,  * 1  ^y.  v3  Yj+T'pl 


(3.53) 


+ 


|  jV-  A  1A  ^J  '^%  ,  fr*-  *  ?  V 


.-"*  2  "  Vjfc I ("'  +  f  ^V(ar  +*-u-ir  tt 


T  AJ(h^4"^)(d°i,^' 


J   ) 
(3.54) 


+ 


A/  \*/^     r  ?  ,W/    *  /-' 


*)kp'+hM<itA*]yl 


TP+h.hPVti!* 


O  o  O 


(3.55) 


h  t'.iW   *.,01\     r.W 


(3.56) 
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-  _  3  — 

+    Tl    (l+fj)    '    J    Qr(;+f3   )    ,  (3.57) 


o 


(3.58) 


p 


itkfa*+w+'K*]K<^iwx>+ti)}- 


(3.59) 


Now,  let  us  establish  the  displacement  boundary  conditions 
which  are  obtained  from  Eq.  (3.44),  i.e., 


o^V 


(  T-  -  V-L  )  &  6 L   J  A    =  0 


H2 


But,   from  Eq.    (3.45),    the  expanded   form   is 
->A  A 

A  (\rrvi)  &i  JA  +\   (vr  nrL)  6S,  JA 

+  1    (fy  -Viji&l   dA  =  0  <3-60) 

for  the  three  layers.   By  again  applying  the  fundamental  theorem  of 
calculus  of  variations,  the  displacement  boundary  conditions  become: 


on  the  top  surface, 

'V-   -     VL        =  0    J  (3.61) 

on  the  bottom  surface, 

V"l    -  1T;      -      0   j  "         0.62) 

and  on  the  edge  surface, 


If-     -     VL      =  0  ,  (3.63) 


ll  II  ~s 

VL  -   Vi   =   0. 


U3 


3.3.3.      Strain-displacement   equations 

The  strain-displacement  relation,   Eq.    (2.7),    is    a  direct  result 
of    the    fourth   integral   of   Eq.    (2.2);    therefore,    consider   the  equation 


[*ij  -{{nlj*  vj^  *  v%  irr/y)]Ssij  d\f  --  0. 

oV  (3.64) 

The  notions  employed  in  Sect.  3.3.1  are  also  applicable  here.  First, 

Eq.  (3.64)  is  written  in  terms  of  the  three  layers  of  the  sandwich 

plate,    i.e., 


vK-2-(VV''r7<M^VK 


o 


+J-  [  *ij  ~    iifilj    *    Vjk    *    VrlL    Vrlj)]llUJ\f        (3.65, 


y 


Lj  £.      X         L,j  J-C  L  '    J 


Although  the  integrals  can  be  combined,  there  is  no  advantage  to  do 
so  in  this  case  because   8'S  -1,   KSL:),   8"S  ^     are  independent 
quantities.   Therefore,  let  us  consider  just  the  first  integral  of 
Eq.  (3.65)  and  rewrite  it  by  splitting  up  the  tensor  as  before; 
then,  we  obtain 
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/,  ft  ^ " i  ( 'V<U +  V*  +  'v %  'V*L  *  'Vj,<  irM)] b's"f 


.°*P  2  <  u*ip   "fiU    u  k  "tip 

0 


J3 


+ft*r  i  (z  K*  +  b\s  K*  +  Ki  '%3ps33}  Jt 

(3.66) 
The  second  and  third  integrals  of  Eq.  (3.65)  will  be  .the  same  except 

the  prime  will  be  replaced  by  a  bar  and  double  prime,  respectively. 

Now,  the  fundamental  theorem  of  the  calculus  of  variations 

can  be  applied,  and  the  strain-disDlacement  relations  for  the  top 

layer   are 


\ =  {{'%Y  V-< +  ^%  ty +  K*  K(*)> 
'K<3  ~-  2  [K* +  K* +  'v*L  Ks  +  Km  Ks), 


(3.67) 


k5 


Note  that 

therefore,        'y  ,  =  'v„      as  expected. 

Too  Y3cf 

By  substituting   the   displacement    functions    Eqs .    (3.1U)    into 

Eqs.    (3.67),    the   strain-displacement   equations    can   be   written    in 

the   form 


where 


/  /  I  -  3  i 

^33     =     X3 

+($f(j%-ir%)(rrirfa/fi) 
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/     r   / 


X/3  -  ilr<ip  +  ''ffiU  * d  %  W//s *  *  4  r7  % 


a  /  -.//  '  */  ^  <. 


*  tf  j„  faf  *  Ujy  fs9<  *  f  (ip%-rl<  )  fy/f 


i  ^r 


^.7,    ,      '      ,    i  ',.87     -x   A 


(3.70) 


:^    "   2   ^  %    W/^   *   ^W   K/3  J 


(3.71) 


A/3  2 


♦I 


2 


^  *  ^3,^    *   ^  7^    ^     *     tl 


<  fs 


?*>*-'?*,*) +  i(f*L-r%)h 


fsW   '    tys,*)    ^3 


(3.72) 


&3  =  z  (  Pr  f  U  f  ^  K°<  *  ^w  J  , 


(3.73) 


1,1  I       V  I 


hi   ~   1   K2  ^3   +   f     f     +    f3^3   ) 


(3.74) 
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The   strain-displacement   equations    for   the  middle   and    lower 
layers    can   be   obtained    in  a   similar  manner.      However,    because   the   equa- 
tions   are   rather    lengthy,    it    is    sufficient   to   observe   that   the  equation 
for      "v.  .      can  be   obtained   from   the   above   by  substituting   a   double 
prime   for   each   prime   and   a      -   h      for   each      +   h;    similarly,    the   equations 
for     v.  .      are   obtained    from   the   above   expression   for      'v.  .      by  substi- 
tuting   a  dash   for   each   prime   and   by  setting      h   =    0. 


3.3.4.      Stress-strain  equations 

As  before,  the  stress-strain  equation  Eq.  (2.8)  is  a  direct 
consequence  of  the  fifth  integral  of  Eq.  (2.2);  therefore,  consider 
the   equation 


SiJ-U£+    H)\bt.M    =    0.  (3-75) 


2  v  6*i\     3  Yji   JJ      <l 


In  precisely  the  same  manner  as  before,  the  integral  can  be  written 
in  terms  of  the  three  layers;  and,  because   *'v.  •»   *v-  •>   6*v.  .   are 

J  TL]        1 J         1-  1 

independent  quantities,  the  fundamental  theorem  of  calculus  of  varia- 
tion yields  the  equations 


5  J 


I*      1L 


2-  K  }%■    yrjc 


S  J  =    9   (  Tzr    +     —  (3.76) 


lJ      V"J' 


1  K  iViy       i%- 


k8 


for  the   upper,   middle,   and    lower    layers,   respectively.      Now,    in   this 
work,    stress-strain  equations    in  the   plastic   range   are   desired.      In 
addition,   the  results  must  be   consistent  with  the   theory  of   plasticity 
for  small   strains.     The   following  method   of   obtaining  the   plasticity 
equations   for  finite  strain   is  basically  due  to  Novozhilov   [6]. 

In  the  modified  Hell inger-Reissner  theorem,    it  was   assumed   that 
the  strain  energy  density  is   a  function  of  Green's  strain  tensor,    i.e., 

If    it    is    further   assumed   that   the  material    is    isotropic  as    in  most 
theories   of   plasticity,   then  the   strain  energy  density  can  be  written 
in   terms   of   the    invariants   of   the   strain  tensor 


Z     -   Z    (f,,I2,fs) 


where 
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Then,  the  first  of  Eqs .  (3.76)  can  be  written  as 


where  the  primes  have  been  dropped  for  convenience.   Next,  from 
Eqs.  (3.77), 

hi  </ 

ill    .   i-;  r'J  -I,  f  - 1,  q'J . 

Substituting   Eqs.    (3.79)    into   Eq.    (3.78)   gives 


(3.80) 


J  A 


Wr'-zi^-i^v) 


Similar  to  classical  theory,  a  deviatoric  stress  tensor  can 
be  defined  in  terms  of  the  Kirchhoff  stress  tensor  as  follows: 


or 
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where 

Also,    in  most   theories   of   plasticity  a  yield   criterion   is   assumed. 
We   assume   a  yield  criterion  of   the   form 

T)   Tj.    =    ZK  (3.8D 

where  K  is  a  material  parameter  depending  upon  the  amount  of  strain; 
note  that  Eq.  (3.81)  reduces  to  von  Mises'  yield  criterion  for  small 
strain  theory. 

When  the  yield  criterion  Eq.  (3.81)  is  rewritten  in  terms  of 
the  Kirchhoff  stress  tensor,  we  obtain 

Similar  to  the  strain  invariants,  there  are  also  stress  invariants, 
the  first  two  of  which  are 


4*  J* . 


'     In  terms  of  the  stress  invariants,  the  yield  criterion  Eq.  (3.82) 
becomes 


(3.83) 


2  4  +  I  J,     =   2K.  <3-8i° 
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Also,  note  that  through  Eq.  (3.80),  the  yield  criterion  can  be  written 
in  terms  of  the  strain  energy  density  and  the  strain  invariants;  thus, 


lj 


ti^z 


(3.85) 

In   addition,   substitution  of   Eq.    (3.80)    into  Eq.    (3.83)   gives 

sz*   ,T  JZ*         JZ* 

J''-3ri,'2I'JJ2'^ri3  ■ 

Now,  in  the  deformation  theory  of  plasticity,  it  is  assumed  that  the 
strain  tensor  is  proportional  to  the  stress  tensor,  see  Kachanov  [20]; 
therefore,  expressions  of  strain  invariants  are  proportional  to  corre- 
sponding expressions  of  stress  invariants.   Eqs.  (3.85)  and  (3.86) 
can  be  made  consistent  with  this  proportionality  requirement  if 


~  * 


-3 
and  thus  become 


'-   =  o 

i  I,         '  °-87) 


<JZ*\Z,       ,    2    rZ 


J   -  3  *Z*     'zr  H' 
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Solving   Eq.    (3.88)    for      9T  /?il2      yields 


JZ* 

h 

d2 

-r 

1 

^2 

nz 

\/2Iz 

f 

2. 
J 

rr 

or 

SI2 

=  2 
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K 

where 

(3.90) 


<  -  4  (j  *.'  yJ  -  r*  ri  f  . 


The  numerical  factors  in  Eqs .  (3.91)  have  been  chosen  so  that  in 

classical  theory  S    and  v   reduce  to  the  octahedral  shear  stress 

J        o        o 

and  strain,  respectively.   Now,  Eq.  (3.89)  can  be  written  as 


IZ*      .  ~  -S 


(3.91) 


J'  =  3u,-  **>  i 


or 


* 


Thus,  the  expressions  for   ^T*/5I.   have  been  determined,  and  the 
stress-strain  equations  can  now  be  established  by  two  different  pro-  , 
cedures . 
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For  the  first  procedure,  Eqs.  (3.87),  (3.90)  and  (3.92)  are 
substituted  into  Eq.  (3.85)  to  obtain 


slJ-jU++t,%)3iJ+*%(rif-j;iiJ), 


t 

or 


t  3  £         J 


(3.93) 


o 


As  already  noted,   J   is  proportional  to   I1  ,  i.e., 


J,    ==  3k L 


(3.94) 


where   k   is  a  physical  constant  of  the  material  and  corresponds  to 

the  bulk  modulus  in  classical  theory.   Then,  substitution  of  Eqs.  (3.77) 

and  (3.94)  into  Eq.  (3.93)  gives 

su  =  2  s  //+ 1 ,3k  _  2  5  ) ft  ij  _  (3>95) 

In   addition,    the    inverse   relation   can  be   obtained    from   Eq.    (3.95)   by 
using  Eqs.    (3.77)    and    (3.94);    thus,    the   exDression   for   the   strain 
tensor    in   terms    of   the   stress    tensor    is 

Eq.    (3.93)    is  Hencky's   stress-strain  equation  for   finite   theory  as 
given  by  Gleyzal    [21]. 


5U 


For  the  second  procedure,  it  is  desirable  to  have  an  explicit 
expression  for  T.   (yf  •)   in  order  to  apply  the  equations  resulting  from 
the  variational  procedure,  i.e.,  Eqs.  (3.76).  Consider  the  expressions 
for  d£*/dl.,  i.e.,  Eqs.  (3.87),  (3.90)  and  (3.92),  which  are,  respectively, 

£*-   0, 
6U  K  ' 


From  Eq.    (3.9U),    the    last   equation  can   be  written   as 


hi    =  1  (3k  *  f  %  )  I,  .  <>•»» 

Then   Integration  of   Eqs.    (3.87),    (3.90)  and    (3.97)   gives,   respectively, 

Together,   the  above  equations    imply  that 

2 
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or,  from  Eqs.  (3.77), 

which  is  the  desired  result.   Naturally,  E.q.  (3.98)  can  be  put  in  many 
different  forms,  such  as 


*_    Sa      kr  hv     v/  //    ,         S.  x    kr  h  v      v 

o  v  o 

(3.99) 
or 


(3. 
Observe   that   from   Bq.    (3.99), 


^  =  ^  ah -J*  ( t,  si  si  +  rt  «;V 

W;  _  y     J       J        ^      U       k      S  AS       r      A 


iw-'fhYO!*'!  %+$.*:**), 


-f  T 

o 


or 


-  =  -«*+*'>  iw-^Xd'Whh1'- 


ir«    < 


Because 
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the   equation 


,  ,  ii        JI 


*v-i(- 


zM*ij     *  fc 


becomes 


which  is  precisely  Eq.  (3.95). 

Later,    it  will   be   desirable   to  have   the   stress-strain  equations 
in   the   forms 

ij  tiki      , 


which  will  now  be  developed. 

Because  of  the  symmetry  of  y   , 


Y 


ij         1/      Ik    j  I         ii    jk 


I    (       iK     J*  <■*     J*    \    s/ 

1(9   3   -3  3    '  %t  > 


and 


r*  "  z(3    9     +d    3    hik  V 


Therefore,  Eq.  (3.95)  can  be  written  as 
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S     -  I  (3    3    +  3    3    )  hi 

I   /  o  y   „  vS,  \  /  b    ks        h    kr\  LJ    y 


or 


CI 

Similarly,  because  of  the  symmetry  of   S   ,  Eq.  (3.96)  can  be  written  as 


y'J       1  I'  /     ik    J*  it  jk  \  c 


+  Lf-.  -4  $)0'JnuSii. 


3  \3k       Z   6.  )  3 

Thus,    the   stress-strain  equations   become 

d    '  L         hi  7 


(3.101) 


where 


(3.103) 
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JJ'M       I    l  ,    Ik  j*       U   jk  ,  /  f  J_   _  1    t  )     ij  kl 


(3.104) 

Note  the  symmetry  relations 

V^i  LjJk  Jiirl  klij 

E         -  E         =   £.  -    E      J  ; 

(3.105) 


F    =  f 


Because   of    the   preceding  work,    it   is   desirable   to  have   the 
stress-strain  equations    in  terms   o£  subtensors.      For  examole, 

33  ^33fh      .  ^33f3  ^  3333      . 

S     -  E        ^  *  2E        fn+£        Y33  . 

But  because   of   Eqs.    (3.103), 


«fl<3  *3fh  «333  33/3 

£         =    E  -    E         =    E  =    0.        (3-l06) 

Therefore, 


^3=  2£-""r„, 


(3.107a) 


59 


^33  33  n  3333      , 


where,    from  Eq.    (3.103), 

«fH         S0    /    Hi   pi  o(i 


(3.107b) 


E  t 


0 


W*  fY>-i  (*-'*  W  3*  > 


(]„ 

"  1  3     • 


(3.108) 


3    Y.    +   k   ■ 

From  the  above,  the  subtensor  form  for  Eq.  (3.102)  is  obvious. 

In  this  work,  it  will  be  assumed  that  the  above  stress-strain 
equations  apply  to  all  three  layers  of  the  sandwich  plate;  however, 
it  is  not  assumed  that  the  material  of  each  layer  is  the  same.   Natur- 
ally, as  indicated  by  Eqs.  (3.76),  different  stress-strain  equations 
could  be  apDlied  to  each  layer  which,  indeed,  would  be  necessary  if, 
for  example,  the  middle  layer  is  assumed  orthotrooic  as  is  often  the 
case. 
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2.U.      Buckling  Problem 

In  thi9  work,  the  bifurcation  theory  will  be  used  to  formulate 

the  buckling  problem.  The  statement  of  the  theorem,  quoting  from 

Novozhilov  [6],  is 

The  moment  of  appearance  of  a  possible  bifurcation  in  the  solu- 
tion corresponds  to  the  critical  load.   Hence,  two  positions  of 
equilibrium  corresponding  to  an  infinitesimal  increment  in  the 
critical  load  differ  from  one  another  by  an  infinitesimal 
amount. 

Note  that  the  buckling  criteria  is  associated  with  two  equilibrium 

positions;  therefore,  all  of  the  preceding  equations  are  aoplicable 

to  both  positions  providing,  of  course,  that  the  inertia  terms   f 

are  omitted.  To  distinguish  between  the  two  positions,  a  dot  will  be 

used.   For  example,  the  displacements  corresponding  to  the  position 

which  becomes  unstable  will  be  denoted  by  v.,  and  the  other  (stable) 

L 

position  will  be  denoted  by  v..   Following  the  procedure  of 
Novozhilov  [6],  these  displacements  can  be  functionally  related  as 
foil ows : 


VL      =  V}    +     h '  V-  (3.109) 

where   6V.   is  the  infinitesimal  change  that  occurs.   Furthermore,  we 

will  assume  that  the  functions   V.   are  finite  and  that   6   is  an  infin- 

l 

itesimal    quantity  which   is    independent  of   the   coordinates.      Naturally, 
it   is    assumed   that  similar   relations   apply  to  all    quantities;    for 
example,    the   stresses   corresponding   to  the   two  equilibrium   positions 
would   be  related  by  the  equation 
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rj*  sj+  wj. 


(3.110) 


Corresponding   to  Eqs .    (3.9),    the   displacements    for   the   sandwich 
plate    in  the  stable   equilibrium   position  would   be 


-      -      -3   JL 

V-  -  ull  -t  9   yi  ?  (3.11D 


»„3   „  • 


V-    -    Hi    *    9    'if/; 


By  applying  the  interface  continuity  conditions  as  before,  and  by 
introducing  the  notations 


(3.112) 


ill  =  a-  *   &U£  , 


we  obtain  the  following  exDressions: 

'vL  =  i/i+  lift  -fr)+  e3,fi; 

V;    =    U-L    -f    9    $L     ;  (3.H3) 

which  have  the  same  form  as  Eqs.  (3.14). 
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Now,  consider  the  stress  resultants  as  defined  by  Eqs.  (3.28)  - 
(3.34).  For  the  stable  equilibrium  position,  similar  expressions  hold; 
for  example, 

a 

h 

By  applying  Eq.  (3.110),  an  alternate  form  is 


a 

-3    .  -3 

M~ 

Jb 

which   from   Eq.    (3.28)    is 


a 


Jk 


M 

Now,  since  the  relationship  between   ;M   and   'm '   which  is  consist- 
ent with  Eqs.  (3.109)  and  (3.111)  is 


we   find   that 

Sim**  \h&9SJ6l 

as  might  be  expected.  The  variation  of  the  other  stress  resultants 
are  calculated  in  the  same  identical  manner;  the  results  are  as 
foil ows : 
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S{ 


/of/3  -<<p        ,1      c<0  .  , 

w  r  ,  m  \  m  r  ;         (  r6 

,     «rtfl       _  *&      ,,    o/£>  I  I   /  , 


^u^fAl  AdS-3., 


j  \ 


iif 


(3.114) 


J  ^   *  C  ^       f    M        J     (3.115) 

Sim33,  *", '**}•£{& J  ^faj  Jff3. 

h  C  d  (3.116) 


6 


',./'    ,-'' ■',.,'}-{•.: J  '..J    -J    -J    «'i _  -^ 


'bL-h 


W  }  w    }  l/v  ur  -  ur  p  ur  -  ur    7  nr 


i   ' 


> 


^ 


i* 


v. 


e    ,<?  e 

e       #  e 


/j 


(3.117) 


ft 


-3 

9 

//  3 

8 


nn 


Z 

h 


i 


'p',t'(    ] 


\'Bl+yL-;p\  v\c' 

s'       «•'      S'         -i       ~. 


b   +  r 


+  2C     -/       ) 


v  z 

(3.118) 


r,  v./-^ 


j 


(3.119) 


The    quantities      6e  are   calculated  by  substituting     Eqs.    (3.110) 

and    (3.112)    into  Eq.    (3.25);    thus, 


6U 


£'3  -  (<s"+ s^H&i-m'l,*  07% 


+  l 


i(?ir+ i?1/, -?'/,  -  i'r 


+  *('r%+si?%)]+(liatu.u)(s' 


+ 


V ' *  *'$'), 


t        -    £       ■+   <U     [  b,  f  u   L  -t  | 


fu-.ir  +  zCr  It-} 'It 

+eYk}+h22(bLST'f)  +  's"lt>dll< 

T    higher   order   terras, 

'fl3  =   'f13      y  i3 

LL  E-       ■+      0  G.         j-     higher  order  terms. 


Therefore, 


(3.120) 
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and   similarly 

$ec3  =  h!i3[hcY  +H  If-j  (ft%  -  Yh 

Z3  "    LI    1       c  u33  /  (l     "  i  \        ii  Y3  r   .  -  i 


h  /  (  -  <■!        c  *  i  i    \      ~  3  r  "_  i  i    7      L  H  ,  "-  . 

(3.122) 


3.U.I.      Field  equations 


With  the  inertia  terms  fL  omitted,  Eqs.  (3.35)  -  (3.U2) 
represent  the  equilibrium  equations  for  the  unstable  position.  By 
replacing   all   of   the  un-dotted   quantities      ('N      ,      u   ,    etc.)     by 

0" 

dotted   quantities      ('n      ,    u    ,   etc.),    the  equations   then   correspond   to 
the  stable  equilibrium  position.      Thus,   by  using  the   form  of   Eq.    (3.109), 
the   first  equilibrium  equation,    i.e.,    Eq.    (3.35),    for   the   stable    posi- 
tion can  be  written  as 
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+  p   +  &  P   +  p   ■+  $  P' 


Since  Eq.  (3.35)  is  still  a  valid  equation  that  applies  to  the  posi- 
tion that  becomes  unstable,  it  can  be  subtracted  from  the  above  equa- 
tion. Also,  higher  order  terms  such  as 

are  neglected.  Therefore,  Eq.  (3.123)  reduces  to 
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f(V*+ 1  '/')b$%  *  (it***  hf)  Y, 


r 

* 


i- q's'y  * sy y}j  - sy*  w* s>'=  o. 

This  same  procedure  is  applied  to  the  remaining  equilibrium  equations, 
i.e.,  Eqs.  (3.36)  -  (3.42);  then  by  cancelling  the   6,  the  buckling 
equations  are  obtained  as  follows: 
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•tft*a*W)uir+ (!»»** 


"+•«")($+ n) 


^{W-V'hM^l^ll^'^)^} 


n 


*  I'm"- \  'h")  'n * (*"-  { y> r% *  fVi  |  n- 


+ 


-8    -<* 


V        —  V      /'of- 


ff    *yf"}l^'P^r*?'     =    0; 


(3.12U) 


Xk^-'^)^]uv[i(^-^h^i(^!iv 


-t 


/>  l/',^  « 


'r+V')+R*]FlA(if(!»"+tf')+k>"]}''l 


t 


4V'-f>';H  +\{'^-\^)Y'h 


if 


~W4^)Yh  -  U^'+iWr'/r 


-t 


h  /  '  6  '  o(      >>  $  »_ 


f> 


QU" 


t 


3S,    -Y 


-fS^r  +"■%)  *  (f%-  H33)p"+(i%-*')y 


-t 


h     /   VX       "->*' 


P  -  P  -)+C 


-4 


o 


(3.125) 
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W-  \  W)(0</t  <  |  p\t )  f  (^_  tyr)(  i 


+  H% +!??%)+ [(I) 'N?U'M^kn'§i 


J 


/ 


+[(0  &-  a***  y>]  fir  -{(w+p  pn 


/ 


'J  f  /7«7       h  /  7  41      W 


-Q[n+un-rir)i-y[t>t  **% 


+  t.  /  rr.^l       ?,V 


(?%-'<?%)! +  (r%-b's)!f>-< 


+  ('t%-W3 


)'^    -     |    '?+    CJ     -    0    ;  C3.U6) 


♦/Yl/v*  a  v,  v'JM+lw+yz)}) 


p 


Y 


(i%-*")'f'+±y+c'  -  o-    (,.,,7, 
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'n?*  n"*  tf>  )Q„i 


3rf 


/  '   ■ 


f  ay    -fit 


+  //>'    +  /n 


<")**, 


r 


*l I  (  W-  W)  <  M <»]  $,}Y    +( *,#.  I  tf<) > 


J;/ 


+  [U«fir-"/")+*/",]ps,r   +('-"-  *  '->'" 


+  I  '    f      H    '    Pi  \  ' , 


r 


Lfi*u  h  : 


3>/ 


+  W+ §**")&,,  -  (V^V^ 


+  &&  +Qfi$3  +  '(fh  +  'fo+n)+f(i+%) 


+ 


"  ,6  /         /' 


ro+hvlfi  *  'p3+p3  +  y  -  o-  „.MJ 


{[k^-'^)*^]uhf*[i(^.u,^ 


y  /  - 


W         Li 


'J,f 


+ 


f 


kf/i.A* 


wy+W)+k<"]$m  ^(Mf-ivy* 


-k(Mp+kliP\ 


+ 


ZKM<   +  $*<■')?„    -  1(^,1  ^<)y 


2Wk-Yk+f'o*h)-yo+h)h 


%/ 


i>r" 


-  ^ 


(3.129) 
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W-WHOs^l^rm 


nMtf-tW 


+ 


+ 


+  f*('+fa)] 


K 


L  -  QY[tJ3,r+i(?hy-psJ 


+  (i%-«")('*h)-?    'p2+'e3    -    0S    C3.no, 


*(i%-^X'*h)  +  i>'*'^  =0-  (3-13I> 
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3.^.2.  Boundary  conditions 

Here,  the  general  procedure  employed  is  the  same  as  that  used  in 
the  previous  section.  All  the  boundary  condition  equations  of  Sect.  3.3.2 
are  applicable  to  both  equilibrium  positions. 

First,  consider  the  stress  boundary  conditions  on  the  top  sur- 
face which  for  the  unstable  position  are  given  by  Eqs.  (3.51).  For  the 
stable  position,  Eqs.  (3.51)  become 


2P*      ~    zPL        =      °  °'132) 


'cc    -  'cc    =   o 


From   Eqs.    (3.i;9)    and    (3.110), 


or 


1        L 

iP* 

P 

I  '  L3             i      L3          e  1    1 3 

=    E      =   E     +  oz     ; 

• 

If* 

=  'pi  +  *  K , 

a3    a    / 

AT        6     -   -r  +    h 


AT 


e-k  +  'h 


(3.133) 


iP'  =  2P'         :r"' > 


+    if 
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where 

hzP \      =  S  Z«     ,       AT         6     =  J   +  ^  , 

By  substituting  Eqs.  (3.133)  into  Eq.  (3.132)  and  then  subtracting 
Eq.  (3.51),  we  obtain 


For  convenience,  the  6  can  be  cancelled;  the  stress  boundary  condi- 
tions for  the  buckling  problem  become 

for  the  top  surface, 


for  the  bottom  surface, 

ii      L  a     I  h      i  ii     i 

and  for  the  edge  surface  (neglecting  the  higher  order  terms), 


(3.134) 


(3.135) 


7U 


I     ci  -°<  W 


O  O 


-^{(.'Ni^N^V^U'l, 


+  ik,t,+Mf*+h?'xst1+u.%)  *  (Vr-/V)H 


f 


-  (W  f  v;  ^7,  *</£<*  <?/*#  **  7^ 


■< 


+YY+i'r+yv)> 


(3.136) 


o  o 


h?-j^{{\W'-'Ni«)+Mf')i?t 


r 


-t 


h/tfif  i  fit 


2  (  M    -  /ri 


)  +  **](&  i%)+[(if(P 


+w*)+^]p%  +[(ifa"+^)*i"]  rh 


2{/Vl    ^  j/V     J?  It         2(W     +  J  M1    J  ^  » 


/ 


(3.137) 
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[(z)    '^  '  h/M^    +   k      ]  '^7r    J  ,  (3.138) 


r 


(3.139) 


$'  $s+f'o*h)+ t'o+fry+yi'+hl 


(3.1U0) 
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1&-.W*.*,  {iiw-  w'htns*, 


r 


f 


uiii'AW).R"]i„t  *  las-  iw&t 


^ifa"**")**"]?*.*  - 1  (W+iwk. 


h  fit    fit     h  " 


(3.141) 


V- 1  'J,  *  It'uP.hW 


t-i*-M(M 


h     - 


*P)(Vs,r  +  I?„r) 


(3.142) 


i3  ,    h     »  3 


\v-jywr*r*f'Khr-i**) 


+(»"+$  *"Xe*r'iM  +  &%*§  tfh 


2 


(3.143) 


77 


where 


I     '  L  I        L 

P     c 

0)     II  i  II,,  L 

P    c 


>=H 


'  ^  L 


-3        •  ^ 

*    T     3     ^* 


II     (. 

— ' 


*    ) 


-3 


.} 


I     L 


_  t         ii    L 

o  >   o  I  o 

S\.i    -i  .«■ 


•\ 


1  ~s  L 


-3 
AT   8   A 


C 


'~h     'l 
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In   the  equations,    the   quantities      6s+      are   the  variations   of   the   pre- 
scribed stress   vectors   that  occur  during  buckling. 

Next  consider  the  displacement  boundary  conditions  which  for 
the  unstable  equilibrium  position  are  given  by  Eqs.  (3.61)  -  (3.63). 
For   the   stable   position,    Eq.    (3.61)    becomes 


Vi 


v-   =  o, 


or 


v-tiVi)-  ('%  +  *%)*  o. 


(3.146) 


Subtraction   of    Eq.    (3.61)    from   Eq.    (3.146)    yields 


S    V;      ~     S    V;        =      0 


/ 
L 


After  cancelling  the   6's,  the  displacement  boundary  conditions  for 
the  buckling  problem  become: 
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on  the  top  surface, 


/  i  ~j 


VL    -     VL      "     0  -}  "  (3.1U7) 

on   the   bottom  surface, 


U  //~J 


VL    -     V-      =      0\  (3.148) 


and  on  the  edge  surface, 

'K-  -  'vi  -  o, 


V.     -      V.        -      0  (3.1U9) 

L  L  1 


Vi  -  v,  -  o. 


In  the  equations,  the  quantities   8V.   are  the  variations  of  the  pre- 
scribed displacements  that  occur  during  buckling. 


3.U3.   Variation  of  strain 

Naturally,  when  the  displacements  undergo  the  incremental 
changes   6V. ,  the  strain  tensor  also  experiences  an  incremental  change. 
To  calculate  the  incremental  change  of  the  strain  tensor,  first  observe 
that  Eqs.  (3.67)  apply  to  the  unstable  position,  and  the  same  equations 
apply  to  the  stable  position  when  all  the  quantities  are  dotted.   For 
example,  in  the  stable  position,  the  first  of  Eqs.  (3.67)  becomes 
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K/s  =  i  (ty  +  V/^  %  'Vrip  +  ^  ^3 >/?  ) ■        (3 .150) 

Substituting   Eq.    (3.109)    into   Eq.    (3.150)    gives 
or,    by  using   Eq.    (3.67), 


where 


Naturally,  the  6  can  be  cancelled.  Also,  through  Eqs .  (3.14)  and 
(3.113),  r^  can  be  written  in  terms  of  the  displacement  functions 
Ucr'  ^a'    *a'   etc*   In  this  manner»  ^  obtain  the  following  equations: 
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h?  -  f<f  *  '%  -  ('Vfy 


■) 


K  -  K  -  K  *  83;r<3 ;  o.»» 


Z7    =  T7    • 


where 


0/  «V9       '2 


(3.152) 
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,r<P  =  i  { pjp  *  iPu f  YL I  otif  *  -z  (f,i  -  f,/. )] 


*WL-  YlMi/ill^  +\iwf-  hp] 


+,uY//h'^   " 


(3.153) 


J 


'n  -  I /'.'/i  ' 


\f  i  Crt  hf  *  f  I  hit  *  'f^hf  *  h<  K/3 1. 


(3.154) 


/    <    '- 


(3.155) 


f<3  -  2  (K*  +fK   ft  t  ft  ft  '  $3>4  f3  *  fs,S  fs  ) 


5 


As =  jp 


/    * 


3  +    Y '  <Py  +  p3   <p3    • 


(3.156) 


(3.157) 


82 


The   above  equations   can  also  be   aoplied  to  the    lower   layer  by  replacing 
the   primes   by  double   primes   and  the      +  h/2's      by     -   h/2's.      Similarly, 
for  the  middle    layer,    the   primes   are  to  be   replaced   by  bars   and   the 
h/2's      are   set  equal    to  zero. 

3.U.U.   Variation  of  stress 

To  calculate  the  variation  of  the  stress  tensor   6s    due  to 
the  variation  of  the  strain  tensor   6F.  .,  one  can  start  with  Eq.  (3.101), 

L  ] 

i.e. , 

the   incremental   operation  on  this   equation   yields 

SxLJ   =    tLJ  hrkj[    +    ^        S£CJ  .  (3.158) 

It    follows    from  Eq.    (3.103)    that 

ijkl         ik  jl      a   jk    2     ij    kl  \</S.) 


Because 


/'    0  Jo  J0     \      b  X> 


*K      K  /   K 


and   from   Eqs.    (3.91) 
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0  7  Y     \  S  "*»    S    '       '     r    5 


Eq.  (3.158)  can  be  written  as 

p 

(3.159) 
The  above  is  an  expression  for  the  stress  variation  in  terms  of  the 

strain  variation  and  the  strain  tensor.   However,  for  most  buckling 

problems,  it  is  more  convenient  to  have  the  expression  in  terms  of  the 

stress  tensor  rather  than  the  strain  tensor.  From  Eqs.  (3.102)  and 

(3.104),  it  follows  that 


kl        Yk£, 


iLt      ~     Vv  /       vS 


where 


(3.161) 
By  substituting  Eqs.  (3.160)  and  (3.161)  into  Eq.  (3.159),  the  equa- 
tion for  the  stress  variation  can  be  put  into  the  form 


bJ-  MiJkUr> 


k£  (3.162) 


where 
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i  ikl 
Note  that  M      has  the  symmetry  characteristics 


likl  1,1k  jtk£  kiii 

M        =  M J       -  M  =   M      J 

In   classical   theory,    the   quantities 


(3.163) 


(3.164) 


SSo  7     4 

are  defined  as  the  tangent  shear  modulus  and  secant  shear  modulus, 
respectively,  and  can  be  obtained  from  a  simple  stress-strain  curve, 
see  Teodosiadis,  Langhaar ,  and  Smith  [22]. 

Naturally,  the  stress  variation  tensor  can  be  written  in  terms 
of  subtensors 

&A,     =  M         MY<+2M        6/L+M         6/ j*  (3.165) 


,     33  3in  3313  3333      n 

<U  =  M     8Pn  +  2M     trr3 1  m     iP33 

where,    from   Eqs .    (3.106)    and    (3.163), 
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M 


M        =  E        + 


I  $So         So   \       *3     f3 

35az\  hi      z)6  6 


(3.166) 


M 


4333 


=    Ml 


33  d  3 


2  liS°-$)(s33-s)J3 


tJ333 

M        =  £ 


3333    A.  /hs0     s0 


+ 


1-1)0^/ 


35;  (  br0      I 
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3.5.   Modified  Buckling  Problem 

To  solve  the  buckling  problem  as  it  has  been  formulated,  one 
must  first  solve  the  fundamental  equations  of  Sect.  3.3  for  the  pre- 
buckling  displacement  functions   u.,  1/.,   etc.;  next,  these  functions 
are  substituted  into  the  buckling  equations  of  Sect.  3:4,  which  are 
then  solved.  This,  obviously,  is  a  complicated  and  tedious  procedure 
which  is  rarely  performed  even  for  simplified  theories,  see  Bolotin 
[23].   Instead,  the  assumption  is  usually  made  that  the  influence  of 
the  displacements  prior  to  the  loss  of  stability  has  a  negligible  effect 
on  the  buckling  phenomenon  and,  therefore,  can  be  neglected.  Thus, 
the  buckling  problem  can  be  simplified  by-  omitting  the  displacement 
functions 

in  the  equations  of  Sect.  3.4. 

3.5.  i.  Field  equations 

In  view  of  the  above,  the  field  equations  of  Sect.  3.4.1 
reduce  to  the  following: 


(rft+aV.  tffi)  +  (M  Rf+  'tit*)  0% 


♦<  f  ^->W'jf7/*  C*1*-  IWM 


^{'NPfylltifF+tf  $'+$% 


\ 


lal  ~°<  "no(  n 

+  P      +     P        +     P  =       0    ■  (3.167) 
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t 


nifty*  w+RnH+Uw-  swp, 


\(^jnn+hiC4P-Qr<p')} 


f3 


I     0<      „    j  _  oi 


Q'n * (f%-M33)^-f,  |  {Kpjtf-  o 


•) 


(3.168) 


UV-  k'Jfi*  /ijPAWn/,7<i  +  A  .7.* 


'jMr-t  [M  -  j 


■h\2i 


(n+zn) 


mv^-hM^KnH—zQ'^}) 


33  \  '* 


-M0%+i(n-na+(W-*"y? 


y-lv+V'  o; 


(3.169) 


'•m* +■[&,*+( m*+  i  A/m  <?%  - 1  $  % ) 


+ 


h  \2  iJi    r 


-i 


k-w+wlH  +  i&'P''}!, 


p 


-¥[n-hn-n)h(H-«33)* 


of  ",-,    </ 


¥  +  ir+c\-  o, 


(3.170) 
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+  ('Ai<s'-^nh,-+(^§^r)Kr 


v 


'     A     i 


,- 


It   A     1 


J 


'3         -3        //    3 
-t   P     +   P      +  P         *     Ox 


f 


(3.171) 


(|  (N^-'N^M^IO^+l^-  I  AT)! 


2 


3,r 


-/• 


3,Y      z 


'P 


(3.172) 


-\'m^'m'»)(u,,^^ 


-   b.  I'n?  ,T  '  0 


3,f    '    ifhl)    ~    2(<f?3    *f 


(3.173) 
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*  [{[fiiP*  i  tf* '**]  hr }  ifi  *  (H-'*")h 


(3.174) 


In  these   equations,    the   quantities 


are   still    defined  by  Eqs.    (3.28)    --(3.34)    and    (3.114)    -    (3.119), 

13 

However,  the  quantities   6e    now  reduce  to  the  following: 


ie'-SJS'+Piify+ltsp/r-sM) 


-3 


S     hf%]       +    S33  bfl    i  (3.175) 


L3  ,        L3         -  VI   ,        _/,  -J      - 


he  "  =  U '%  S  13(it/%+f&!plJ,) *  S33h$  ' ; 

(3.176) 

-3      "/'in  "33       a    ; 

f    9       bfLlt]    f      S        6ft        .  (3.177) 
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3.5.2.   Boundary  conditions 

For  this  case,  the  equations  for  the  stress  boundary  conditions 
on  the  upper  and  lower  faces  are  the  same  as  those  given  in  Sec.  3.5.2, 
i.e. , 


on  the  top  surface, 


2<    *    ~  z'         =    0 \  2^ *   ~  2^       ~     ^5         (3.178) 


on  the  bottom  surface, 


",-s  L  "n  l  II     i  h 

l(~*+tC  U,  (3.179) 


However  on  the  edges,  the  stress  boundary  conditions  now  become 


(3.180) 


(3.181) 
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I-' 


-pAhtf-itf+OfLh 


z...  .„..  ,..  ^vm 


+[({f'^-bh^nn+(T?-W)y 


j 


(3.182) 


'^,  h   m      -   A 


-/- 


(3.183) 


3       -3     ,,  3 


p+?+p  =  *fi{(^^y<)uhY  +(k*.kn*fat 


-h 


+ 


[k«fir-^)*^']hr+(^l^r)h 


r 


cfki-Jfr  +  Ph+P+p+y}, 


(3.184) 


Y 


~l(^l^)h,f+i(^-Qh) 


+  11  /  UP-  ".J 


yi+T'h  +  v), 


(3.185) 
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:i3-l^%^{('^-i'^)(0h^ii3l,) 

+  (Tf-?'q')$3}  ,  C3.») 


The  definitions  of  Sects.  3.5.1  and  3.1.2  apply  to  Eqs.  (3.178)  -  (3.187). 

For  the  displacement  boundary  conditions,  the  equation  of 
Sect.  3.4.2,  i.e.,  Eqs.  (3.147)  -  (3.149),  hold  for  this  case  also. 

3.5.3.   Variation  of  strain 

By  reducing  the  equations  of  Sect.  3.4.3,  the  variations  of 
strain  for  the  upper  layer  is  given  by  the  equations, 


^Vvv^v^A^ 


^=<^=-<^^V,3, 


^33    =  <5/733  , 


(3.188) 
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where 


/^J*2  (£*/,■*  V^> 


2'  rfjfl        0  , 


(3.189) 


A     =      ^  -  ' 

As  before,  the  equations  for  the  middle  layer  are  obtained  from  the 
above  by  replacing  the  primes  by  bars  and  by  letting  h/2  =  0,  and  the 
equations  for  the  lower  layer  are  obtained  by  replacing  the  primes 
by  double  primes  and  the  +  h/2's   by  -  h/2's. 

3.5.4.   Variation  of  stress 

The  equations  of  Sect.  3.U.U  remain  unchanged  and,  therefore, 
apply  directly  to  this  case.   Naturally,  the  expressions  for  the  strain 
tensor  variations  are  to  be  taken  from  Sect.  3.5.3. 


4.   SIMPLIFIED  LARGE  DEFLECTION  THEORY 

As  was  mentioned  in  Sect.  2,  different  theories  can  be  developed 

from  the  variational  principle  by  assuming  a  particular  form  for  the 

strain  tensor.   In  the  bending  of  plates,  the  displacement  component 

v_  normal  to  the  plane  of  the  plate  is  usually  much  larger  than  the 

disDlacement  components   v   in  the  plane  of  the  plate.   For  most  mater- 

a 

ials,    in   fact,   the  displacements      v       are   often    infinitesimal    quantities; 

O' 

therefore,  the  nonlinear  terms 

can  justifiably  be  omitted  from  the  strain  tensor  to  obtain 

for  the  coordinate  system  developed  in  Sect.  3.1.  Then  the  variational 
equation,  Eq.  (2.1),  takes  the  form 


y 


jSJ (Vii-+  Vjk  *  tr5ti  vhJ)  dv 


A  r,  vL  j  a  A  sj  ( n  -  vc )  j a] 


O^S  o^V 


9U 
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which  can  be  put  into  the  form 


o     5 


A 


from  which  the  fundamental  equations  are  obtained. 


4.1.  Fundamental  Equations 

Obviously,  the  fundamental  equations  can  be  obtained  from 
Eq.  (4.2)  in  the  same  manner  that  the  general  equations  of  Sect.  3.3 
were  developed  from  Eq.  (2.2).   However,  through  the  results  already 
obtained  and  the  variational  principle,  an  alternate  method  is  avail- 
able, namely,  to  deduce  the  equations  from  those  of  the  general  theory. 
Because  the  alternate  method  is  instructive  and  fairly  simple  for  this 
case,  it  will  be  employed. 

4.1.1.   Equations  of  motion 

From  Eq.  (4.2),  the  equations  of  motion  are  a  result  of  the 
equation 
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Comparison  of  Eqs.  (3.16)  and  (4.3)  shows  that  the  difference  between 
the  equations  is  the  terms 

■        L  (siJv%  )l  Svr  Jv-\(sh"3  v3ii)l  hvr  JV 

(SlJvV;)l.    Sir    JV  (...) 

which  must  be  omitted  from  the  equations  of  Sect.  3.3.1  to  obtain 
the  equations  of  motion.   Note  that  the  terms  to  be  omitted  are  asso- 
ciated with  only  the  6v^  quantities,  not   6v3 ;  therefore,  not  all 
of  the  equations  of  motion  will  be  affected.   Formally,  the  equations 
are  reduced  by  letting  v01     =  0  when  v°     is  multiplied  by  6v^. 
For  example,  Eqs.  (3.18)  and  (3.19)  show  that  the  quantities  A 


a 


a2 
and  B  ~  become 


/■P-j*  8-o.s-" 


33   „_J   „33 


while  the  quantities  A    and  B    remain  unchanged,  i.e., 
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When  the  reduction  method  is  applied  to  Eq.  (3.27),  it  becomes  apparent 
that  the  quantities 

"  >    f     >    f    >    f 

must  be  set  equal  to  zero  in  Bqs.  (3.35)  -  (3.38)  and  Eqs.  (3.39)  - 
(3.42)  remain  unchanged.  Therefore,  the  equations  of  motions  are 
as  follows: 


N   +N     t  N    )\  +  p   +  p    +  p      =    CL  *  CL    f  CL     ;  (4.5) 


{{(N^n^n.-Q^U'p-'p)^* 


[CL   -  CL  )   +   /YH 


_        f)      /  I  ai         a  U    \  —    at 


(  it*!6-  h  a/*P)I  -'n       ~h   '^    'n4         h  '  4      <  << 
"*/P^h   "A,4fl\l       "^     h    »  4      "4       h    "4      "     <y 
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+  <f(/+ifs)-t  Qf(t  +  p3)  f  Q'O+'fj)} 


f 


'3  -3  »     3  1     3  -3  *    3 

+  n   +  o    ■+  b     -    &   +  a    +  a 


p  +p  +  p 


(4.9) 


'h\2/isfit  '..fits     -Ml-       -tLf/flF+h/i/P)! 


fs,/ 


A*  - 


+ 


h    /  /  3        "3 


p  -  p  )  *  £ 


-  3 


—  3 


(4.10) 
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h     i   3        i      3 


_     Jo     a  3        n     3 

-    j    CL    +   nn     •  (4.12) 


and   the      E   "       quantities    become 


t     43  '43  -*3  -43  "43  "Q43 

E       =   j      >        E        -  o       ;        t        -    O       , 

E     =S     {u.3j+8  fht)t$    (/+f3)j  (4-13) 

'i   33     ffc/3f-  A/-  *  >     ^J"         7      oJJ/       "       ^ 
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4.1.2.   Boundary  conditions 

For  the  surfaces  upon  which  the  stresses  are  specified,  the 
boundary  conditions  are  determined  from  the  equation 


l[SJ,-/»lSJr(ii  +  iiviir)]&ViSA  -  0 

■J  A  J 


(4.14) 


JA  J 


which   is   a  result   of   the  second    integral    of   Eq.    (4.2).      Following   the 
procedure   of   the  preceding   section,    the   difference  between  Eqs.    (3.43) 
and    (4.14)    is 

ja   °  J  r        l  j^  o    j 

a     S  /-  o      S 

=)    -  M.  Sjr  Vlr    1%    J  *  i  (4-I5) 

and,  therefore,  these  terms  are  to  be  omitted  from  the  equations  of 
Sect.  3.3.2.  Note  that  the  reduction  procedure  is  identically  the  same 
as  described  in  the  preceding  section.  Then,  with  reference  to 
Eq.  (3.48),  the  equations  that  change  are  Eqs.  (3.52)  -  (3.55);  the 
stress  boundary  conditions  are  as  follows: 

on  the  top  surface, 


y.-p-o,   fr'p-o% 


(4.16) 


on  the  bottom  surface, 


fc+,y  =  o,     ,'cl.  +  ;c  =  o ; 
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on   the  edge   surface, 


(4.22) 
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Also,  the  equations  for  E   ,  Eqs.  (4.13),  which  were  established  in 
the  preceding  section  are  applicable  to  Eqs.  (4.16)  and  (4.17). 

The  displacement  boundary  conditions  are  a  result  of  the  third 
integral  of  Eq.  (4.2);  therefore,  consider  the  equation 


iAv 


V:   -    -irL   )    SSI     J  A    =    0.  C*.26) 


Comparing   the   above  with  Eq.    (3.44)   shows   that   the  equations   are 
identical,   and,    thus,    the   displacement   boundary  conditions   are   the 
same   as    for   the   general   nonlinear   theory,   namely, 

on   the   top  surface, 


VI     -     'ifl     =    0  ,  (*.27) 


on   the   bottom  surface, 


V~i      -      V'l       =      0}  (4.28) 
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on   the   edge   surface, 


V£ '     -    V;     =     0 ,  <u-29> 


V-     ~      V{      =       0  ,  (^-30) 


2£  -    Yi    -    0 .  con 


U.1.3.      Strain-displacement   equations 

Obviously,  from  Eq.  (4.2),  the  terms  to  be  omitted  from  the 
general  theory  are 


that   is,    all   the  nonlinear   terms  with     v       are   to  be  dropped   from 
Eqs.    (3.68)    -    (3.74).      Then,    the   strain-displacement  equations   for 
the  upper   layer  become 

^33      =     o*33    7 
where 
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(4.33) 


(4.34) 

fa?  =  i  fa>.°<  fa*/*  7  ("-35) 


/  /  / 


The   strain-displacement  equations   for   the  middle   and    lower   layers   are 
obtained   from  the   above  equations   by  the   procedure   described    in 
Sect.    3.3.3. 
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4.1.4.      Stress-strain  equations 

The  derivation  of   the  stress-strain  equations   that  was   devel- 
oped  in  Sect.   3.3.4   is   also  applicable   to  this   case;   however,    in  using 
those  equations,   the  expressions    for  the   strains   are   to  be   taken, 
naturally,    from  Sect.  4.1.3,   not  Sect.    3.3.3.      For  example,    the   stre 
'S°®      is   given  by  Eq.    (3.107)   as 


S     =  t  Yn   +    E  '      ^ 


ss 


where 


£F    -  j(3*-z$J3 


i 


I     ,    n'i      'J  'ok     >CJ     \i 


a 

I 


4  =  i(^;^-  's;x?t 


and      ' 


V  &t   y33   are  given  by  Eqs.  (4.32). 
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4.2.   Buckling  Problem 

Naturally,  the  equations  associated  with  the  buckling  problem 
can  be  developed  in  the  same  manner  that  the  equations  for  the  general 
nonlinear  theory  were  developed  in  Sect.  3.4.   However,  in  order  to  be 
consistent  with  the  preceding  section,  the  buckling  equations  will  be 
deduced  from  the  general  nonlinear  theory. 

4.2.1.   Field  equations 

The  field  equations  are  a  direct  consequence  of  the  equilibrium 

equations  (i.e.,  the  equations  of  motion  with  the  inertia  terms  omitted); 

therefore,  the  discussion  of  Sect.  U.l.l  is  also  applicable  to  this  case. 

Thus,  the  first  four  field  equations  are  obtained  from  the  general  equa- 

tions  by  omitting  all  of  the  v   terms.   Now,  from  Eq.  (3.109),  the 

ot  (y 

omission  of  the  v   terms  also  implies  the  omission  of  the  V   terms; 

therefore,  by  direct  inspection  of  Eqs.  (3.124)  -  (3.127),  the  first 

four  field  equations  become 


(m^+ /n  r t  m  ^)l    +  P    +  P    +  P    -  0,  (^-39) 

r 


(4.41) 
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The   other   four   field  equations   are   the   same   as    for   general   nonlinear 
theory  and  are  given  by  Eqs.    (3.128)    -    (3.131). 

In  the   field  equations,   the   stress   resultants   are   defined  by 
Eqs.    (3.28)   -    (3.34)   and   the  variation   of   stress   resultants   by 
Eqs.    (3.114)   -   (3.119).      In   addition,    the   quantities      fie  now 

reduce   to 

Sc"3-    S'l/?9      Se*3=  &*\      S'e"*  M*> 
he*  =  SJl'*[  itM  +  ltfa-far)'*3  far! 
+  bA?(t+f3)  +  'jr3[Wj,y  fz(&<p3,y-ifar) 

+  6    S'p3,r]  +  J**  &{p3  •> 

(4.43) 

he     =-  Sz  3(ii3lf  +  §  f3}/)+  S^JJ(/+fi3) 

+  s*3 (bu3^  +  d3b$M)  +  s33 b§3  , 

(  "  33       c  //  Y3  f  -  h  /  r  ",         \       n 3 "  1 

be     ^bA,     [  Lthi-  j  (p3,r-faf)+0    far] 

Hence,    Eqs.    (4.47)    instead   of   Eqs.    (3.120)   -    (3.122)   are   to  be  used 
in   the  definitions   of      P   ,    Eqs.    (3.119). 
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4.2.2.   Boundary  conditiotiB 

As  with  the  field  equations,  only  the  form  of  Eqs .  (3.136)- 
(3.13  9)  change;  the  other  equations  remain  unchanged. 

Therefore,  the  stress  boundary  conditions  for  the  top  and 
bottom  surfaces  are  given  by  Bqs.  (3.134)  and  (3.135);  while  for  the 
edge  surfaces,  the  equations  are 


0  0  0  "A 


o        A     * 


0  £.     o  o      a      x  £  ' 


o  2.     °  o     a  \  Z. 


and  Eqs.  (3.140)  -  (3.143).  As  before,  Eqs.  (4.43)  which  define  the 
quantities   6e    are  to  be  used  in  Eqs.  (3.119). 

The  displacement  boundary  conditions  remain  unchanged  and  are 
given  by  Eqs.  (3.147)  -  (3.149). 

4.2.3.   Variation  of  strain 

Obviously,  the  change  here  is  to  omit  all  of  the  nonlinear 
terms  involving  either  va     or  V®     from  Eq.  (3.150).  Hence,  the 
variation  of  strain  for  the  upper  layer  is  given  by  the  equations  as 
follows : 
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^ 


3 


-3  .2    ,     ', 


V„,  *'"VV  ^V'iA,  , 


-3 


^,3=  V7**  i >  *As  ■> 


(4.48) 


S/7, 


JU 


<V*3     \ 


where 


./ 


f '- I i ifjf  +  ufu  *  j($4+W  PJf- iff/* ) 

'fas  ![(&,<-  h*  )I  h? + 1  C/¥  -  kfi )] 


'fas  z  ($*-  h,  )][*,*  *  i(  he  -  fa* )]}, 


¥ 


(4.49) 


'[n-^+id'h*- philips,  fi  +  'p3,A  *j,f 


h  -it1. 


7 

(1.50) 

(1.51) 


+ 1(  £3^- faj]  ?&[£,,<+%  (?3><-faw)]}, 


(4.52) 
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f.43     '     2    (  $3,4    +  f3    '<p3,4    +  $3    f*>4  )   , 


(4.53) 


3  (4.54) 


For  the  middle  and  lower  layers,  the  strain  variation  is 
obtained  from  the  above  equations  in  the  usual  manner. 

4.2.4.   Variation  of  stress 


The  equations  for  the  stress  variation  are  taken  directly  from 
Sect.  3.4.4  and  are  not  modified.   Naturally,  the  expressions  for  the 
strain  and  strain  variation  that  appear  in  the  stress  variation  equa- 
tion are  to  be  taken  from  Sects.  4.1.3  and  4.2.3 


4.3.   Modified  Buckling  Problem 

There  are  no  new  concepts  introduced  for  the  modified  buckling 
problem.  The  equations  can  be  deduced  from  the  general  nonlinear  equa- 
tions of  Sect.  3.5  by  the  procedure  of  Sect.  4.2,  or  the  equations  can 
be  obtained  from  those  of  Sect.  4.2  by  the  method  used  in  Sect.  3.5. 
In  either  case,  the  equations  can  be  written  directly  by  inspection. 


Ill 


4.3.1.      Field  equations 

The   field  equations   are   Eqs .    (4.39)   -    (4.42)   of   Sect.   4.2.1 

and   Eqs.    (3.171)   -    (3.174)   of   Sect.    3.5.1.      In  addition,    the   quantities 

a    i3  I. 

oe  now  become 

de    -  bA,    ,       dc     =   dA,     j       be     »  b a,    } 


Se^<U33^/^/^ 


~3 '  hi\+ni 

(4.55) 


4.3.2.  Boundary  conditions 

The  stress  boundary  conditions  are  given  by  Eqs.  (3.178), 
(3.179),  (4.44)  -  (4.47),  and  (3.184)  -  (3.187).  As  before,  the  6ei3 
quantities  are  defined  by  Eqs.  (4.55).  The  displacement  boundary 
conditions  are  Eqs.  (3.147)  -  (3.149). 

4.3.3.  Variation  of  strain 

It  is  interesting  to  note  that  the  variation  of  strain  is  the 
same  as  for  the  general  nonlinear  theory,  and,  therefore,  the  equations 
of  Sect.  3.5.3  apply  directly  to  this  case. 

4.3.4.  Variation  of  stress 

The  stress  variation  equations  of  Sect.  3.4.4  remain  applicable 
to  this  case. 


5.      SIMPLIFIED   URGE   DEFLECTION  THEORY  FOR   SANDWICH 
PLATES  WITH  THIN  FACINGS 

In  general,   sandwich   plates   are   constructed   so  that   the  upper 

and    lower    layers    (i.e.,    the   facings)   are  much   thinner   than  the  middle 

layer.      For   thin   plates,   Kirchhoff s    first   assumption   is   usually 

employed,   which   quoting   from  Todhunter   and   Pearson    [24]   is 

Every  straight    line    in  the   plate  which  was   originally  perpen- 
dicular  to  the   plate-surfaces,    remains   straight   after   the 
strain  and   perpendicular   to  the  surfaces  which  were   originally 
planes   parallel   to  the   plate-surfaces.  6  y 

Thus,    for   thin   plates,    it    is   assumed   that   the   deformation  due   to  the 
shear   stresses    is   negligible,    i.e., 

Also,  in  thin  plate  theory,  it  is  usually  assumed  that  the  distance 
between  points  on  the  straight  lines  which  are  normal  to  the  plate- 
surfaces   remain  constant   during  deformation    ,   or 

'33     =     °>  (5.2) 

Novozhilov  [6]  has  shown  that  with- the  above  assumptions,  i.e., 
Eqs.  (5.1)  and  (5.2),  and  the  simplified  strain  tensor  of  Sect,  k, 
as  given  by  Eq.  (4.1),  the  displacement  functions  u1  and   f1  are 
related  as  follows: 
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Then,  from  Eqs.  (3.12), 

Therefore,  corresponding  to  Eqs.  (3.12),  we  have 

~n  — 

VV2»1  few'  I  &,<  A      ij  =  ^  - 1  ^ , 

and   corresponding  to  Eqs.    (3.13), 

I  '/  "l 

—        III*  f)    I  /j   //  \  _  /       /  „ 


(5.3) 


(5.4) 


'3  > 


(5.5) 


(5.6) 


As  mentioned    in  Sect.    3.2,    the   theory  can  be   developed    in  terms   of   two 
sets   of   displacement   functions,   either 

11 L  ;      fl 

or 

/  // 


11U 

For  consistency  with  the  preceding  work,  the  former  set  of  displacement 
functions  are  chosen.   Substitution  of  Eqs .  (5.4)  into  Eqs.  (3.14) 
yields  the  displacement  expressions  for  the  three  layers  as  follows: 

i 

.3 

VL  =  U-i  +  9    fi  3  (5.7) 


5.1.   Fundamental  Equations 

In  Sect.  4.1,  the  fundamental  equations  were  deduced  from  the 
general  equations  of  Sect.  3.3;  obviously,  a  similar  procedure  could 
be  employed  here  to  reduce  the  equations  of  Sect.  4.1.  to  this  special 
case.  However,  care  must  be  taken  if  the  reduction  is  to  yield  the 
same  identical  equations  that  would  result  from  the  direct  application 
of  the  variational  principle,  rather  than  an  equivalent  set  of  equa- 
tions.  For  example,  Eqs.  (4.7)  and  (4.8)  which  are 


f  II    o(6       L     »  c<($  11/        L     u4       ll*t  £)     II  o(  j,      o( 

{Mr+2Np)l-q  +  ip  +  c  =  fa,  +** 


(5.8) 


are  obviously  unaffected  by  the  change  in  the  form  of  the  displacements 
and,  therefore,  should  be  applicable  to  this  case;  indeed,  these 


115 

equations  are  valid.  However,  the  origin  of  Eqs .  (4.7)  and  (4.8)  lies 

in  Eq.  (3.27)  which  shows  that  these  equations  are  a  direct  consequence 

of  the  arbitrary  functions   6'$   and  6" 4  .   Now,  for  this  case, 

a  o> 

*'f J    =  '  (^3,4    *  |  ip3,*)j 

(5.9) 

bfj  »  -  (s  a 3,4  - 1 5  ^«0 

as  seen  from  Eqs.  (5.4).  Obviously  then,  terms  in  the  variational 
principle  must  be  regrouped  before  applying  the  fundamental  theorem  of 
calculus  of  variation;  the  resulting  equations  are 

(5.10) 

-  C    -i-  C      =  j  ( CL    t  CL    )  -  /w     -t  /yrt     . 

Eqs.  (5.8)  and  (5.10)  are  equivalent. 

As  indicated  above,  the  reduction  process  for  this  case  pre- 
sents no  particular  difficulties;  however,  the  process  is  not  nearly 
as  direct  as  in  Sect.  4.1.   In  fact,  the  beauty  of  the  variational 
principle  becomes  obscured  by  the  reduction  process;  therefore,  the 
equations  for  this  section  will  be  obtained  by  the  direct  application 
of  the  variational  principle. 
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The  modified  Hellinger-Reissner  equation  for  the  simplified 
strain  tensor  has  already  been  developed,  Eq.  (4.2);  therefore,  the 
fundamental   equations. are   a  result  of   the   following  equations: 


V 


J        I  (5.11a) 

,   [K-s>.SJ'r(6lr  +  6i,r)]Sv.  JA=  0;  o.i») 

'A  J 

J     (V--V.)    S<>i    J  A   -    d);  (3.XI6) 

Uhj  '  I  (  VJ:  -  Vj/i  +  V3i£   V     ■  )]  6^  M'O; 

[V         J  (5. lid) 


*  I  <^* 


5.1.1.   Equations  of  motion 

The  procedure  that  will  be  used  here  is  in  general  the  same  as 
that  employed  in  Sect.  3.3.1.   Consider  Eq.  (5.11a),  i.e., 

which  in  terms  of  subtensors  becomes 

Jy       r  '  (5.12) 
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where 


(5.13) 
In  terms  of  the  three  layers,  Eq.  (5.12)  becomes 


~b 


+\  [(S^-  )  iv,  *  {A3\  +  ~- )  Sir,  J        ».»> 


Next,  Eqs.  (5.7)  are  substituted  into  the  above  equation  to  obtain 
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I 


Jh 


[    i  4 A  i  '  ^  3 


■'ft 


*i   [S'l* 


■] 


/°o    (o 

C 

+ 


>"-  f 


't, 


h  r  '  *<& 


Vs 


-£  I  r*%* -]}*?< 


lb 


i'-SiW* 


c< '  fi 


f> 


■f 


■] 


Q.    - 

-  I §3-  - 


4 


U  z 


z 


fi 


i- 


?3^ 
'f)3 


/i  if-m 


%[*%  +  "] 'fib"/,*.-]}  ft, 


if'^v-^Hno^-] 


h    r  'L  3 


U    2 


Lb'''lfi*"]}s$s]j6sJ/i  -  o 


(5.15) 
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where 


(5.16) 


—  33  —°<3  _j  _         _jj> 


~     6      (*-*<-    2   fht)    t    6 


By  using  equations  corresponding  to  Eqs.  (3.26),  and  by  using  the  equa- 
tions for  the  stress  resultants  Eqs.  (3.28)  -  (3.34),  the  fundamental 
theorem  of  calculus  of  variations  can  be  applied  to  Eq.  (5.15)  to 
obtain  the  following  equations  of  motion: 


/\/^  /|7^         A/75)/  '^  -°C  "   o(  I    °C  _^         //    Y 

/V    +  N  l  +Nr  )l  +  p  +  p   *  p    -   ol  t  cl  t  a,    5 

^  (5.17) 


2  (a,  -  a,  J  +  /w      ■  (5.18) 


__/?//«•'       'Wi  —    °t 
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f  fa  ,  i  *3     "  *a.       /ota      »  J/3  7  /  '<*       "J     I  ,  iJ    i/  y  v 

/  |  (Nl-  N()-  M?-  M  I  ]/    -t  Q   *  <?  +±(p-p) 

-  C  -   C      =       |  (tf,    -  <C    )  -  *VI     -   Aff     J  (5.19) 

[j(ti' N?)- mf+mp]\ 'j  i '■  q  -q  +z(p+p) 

-V+   C*    «    l^+'of)-'^*'™     -,      (3.20) 

>  6       -6,  _    .      «fi)i        i  3      .3    ii  3         i  3     _3    //  3 

+  q^QF{i  +  Y3)+Q!}\+p+p+p  -  *>+a.+*.} 

r  (5.2D 


f 


■3>* 


Hf%-N}3)(Hf3) 


+  z('p-p  )*  c    =  z(o.-a,)*im   .(sm 

In   the   above  equations,    the   following  definitions   have  been  used   for 


p        and 

L 

/    V        '  V                    /  VJ 

-3  ,^ 

'  mP    j  £     \         (    $       -i 

(9   5    > 

-M            -°< 

_°<3 

-^  -<? 

< 

mP    ,  *£ 

H 5  ' 

0   5       , 

^< 

1    ^       , 

_  3  „  o< 3 

6   3 

/  3 

9  \ 

-3 

9 


*7        r 


2 

h 


i  JN 


K  2  ' 
(5.23) 
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t\ 


'33        -3  '  33 

fE    ,  OE    s 


<'9\ 


-3  -3  -p33        -3 -S3  -3  ^ 


//    3         H.3 


All 


n    33 
I    E 


-3  0  £ 

9  E 


a   3 

\9  ) 


k 
z 


) 


*P  )  * c   ) 

along  with  Eq.    (3.34). 

It  should  be  mentioned  that  when  the  fundamental  theorem  of 
calculus  of  variations  was  applied  to  Eq.  (5.15),  the  quantities  6u 
and      6u_    (also,      6\|r-        and      6iji,)  were   assumed   to  be   arbitrary   inde- 


(5.24) 


3,  a 


pendent  functions.  Rather  than  use  this  approach,  one  could  use  the 
relation 


j.  (•■■)sam3,<Jv  ']  (■■■ )  hn 

.■>"  nil 


3    o"l*    dA 


'(-■)/.  a,  J? 


(5.25) 


to  combine   terms;    this   has   two  effects,   namely,    to  decrease   the  number 
of   the  equations   of  motion  by   increasing  their   order,   and   to  add   addi- 
tional  terms   to  the  stress   boundary  conditions.      Naturally,    the  equa- 
tions  are  equivalent,    regardless   of  which  method    is   used.      In   fact, 
the  equations   of  motion   that  would   result  by  using  the   above  relation 

(and  the   corresponding  relation  between     6i|r_  and      6if„)   can  be 

-3  ,  Of  •} 

obtained   from  Eqs.    (5.17)    -   (5.21)   by  merely  substituting   Eqs .    (5.19) 
and   (5.20)    into  Eqs.    (5.21)   and    (5.22),   respectively,    through  the 
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fna  A  U^tt 


stress  resultants   'q   and  »Qa.     Thus,  the  alternate  set  of  equations 
of  motion  is 

(  N      -f  N      +  N      J  J     +  p    +  p    ■+  p      =     a    -r  CL    -r  CL 


(5.26) 


[i(h+-'A>*)tii*]lr?+§(p-P')  +  c< 

~    j{Cb    -   <£   )  +  /M        .  (5.27) 


{(k*+*****)iiv+[!i(kfit-k't)*Mn?s, 


r 


h  ,  , 


1  (~ p    +p'    +  CL    ~  CL'  J  +  Cr  +    C    -  /TfV   -  ,W 

'3         -3  a  3  i    3        -3         a  3 

+  p     +p      +  p         =      Cb     +    CL      +    CL      ;  (5.28) 

-  ii  ifr**  w)~'m"+  y% + (?%-«%  & 

(5.29) 
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5.1.2.   Boundary  conditions 

Eqs.  £5. lib)  and  (5.11c)  are  used  to  establish  the  boundary  condi- 
tions.  First)  for  the  stress  boundary  conditions,  consider  Eq.  (5.11b),  i.e., 

By  the  procedure  of  Sect.  3.3.2,  the  equation  corresponding  to  Eq.  (3.48)  is 


l\ 


,t  i  (s:-'f)[h-u.^\h^  -  (f-i)  s&„  -!(§'-  *)if„] 


"  ~*    >*3\r..       f._    ,-3    L,,-  L     -3 


(3+syK-lir<-(*^)***,+l(H)¥*<] 


{[jk  -a  w  {(K  -^  s  *;  {(H  -jy  wh 


CI 


-f 


■4  rc 


(5.30) 
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Then,  the  application  of  the  fundamental  theorem  of  calculus  of  varia- 
tion yields  the  stress  boundary  conditions  as  follows: 

on  the  top  surface, 


on  the  bottom  surface, 
on  the  edge  surface, 


(5.31) 


(5.32) 


6  +  ,5  +  S    =  /nA(  N  r+ N  r+ N  r  ),  (5.33) 

oo  C  \Q         o       '       o      6  L  *-    K  ' 


(5.35) 


(5.36) 


5.38) 
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where   Eqs .    (3.28)   -    (3.30)   and    (3.49)   have  been  applied.     These   stress 
boundary  conditions   are   associated  with  the  equations   of  motion  given 
by  Eqs.    (5.17)    -    (5.22). 

For  the  other  set  of  equations  of  motion,  Eqs.  (5.26)  -  (5.28), 
a  corresponding  set  of  stress  boundary  conditions  must  be  developed,  as 
discussed    in  Sect.    5.11.      By  substituting   the   quantities 


Q    +   q        and  q     -    Q\ 


obtained  from  Eqs.  (5.19)  and  (5.20),  into  Eqs.  (5.37)  and  (5.38),  the 
stress  boundary  conditions  are  Eqs.  (5.31)  -  (5.36)  and 

i  3      -J     "3       _  (,*&    -«#  "  *6\  -        ,  f/?//=<7S  »°<e.x     -°<fl}  _ 

+  C  +  C  -/mF-/Mp  +  Q    (i+  ^3)  j  }  (5-39) 

(5.40) 
The  displacement  boundary  conditions  are  the  same  as  in 
Sect.  3.3.2  and  are  given  by  Eqs.  (3.61)  -  (3.63). 
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5.1.3.      Strain-displacement  equations 

From  Eq.    (5. lid),    the   strain-displacement  relations   are   a   result 
of   the   equation 

,  [Itj-I&tlj*  rjli  *  »5,<  % )]  S6'J^  '  °- 

o 

With  the  basic  premise  of  this  section,  Eqs.  (5.1)  and  (5.2)  for  the 
facings,  the  above  equation  becomes 


7 


0 


when  written  in  terms  of  all  three  layers.   From  the  fundamental  theorem 


of  calculus  of  variations,  after  the  middle  layer  has  been  expressed 
in  subtensor  form,  Eq.  (5.41)  yields 

K/> =  i  ( vjfi  *  'fpU  +  fa,*  hr-3,p)  ; 


^33    =     V3,S    +    {    V}>3   V3>3     , 
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After   substituting  the  expressions   for  the   displacements   Eqs.    (5.7) 
into  the  above,    the   strain-displacement   equations  become 

i  I 

hi,    -  -       -      .        -->.- 


Hi)  [zfs^  *  rM  f3,p)  -  ze  (u.3,^->y3^)], 

(5.13 


*  ">3>J  ftp  +f3>4  &J,f  )+    (0  )     <p3i<  pJifi], 

C5.' 


5.13) 

b-t  ^3,f 

(5. 11) 


l  r  -         -  -    _  -3 ,  - 

(5.15) 
(5.17) 
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5.1. 4.      Stress-strain  equations 

The   stress-strain  equations   are   obtained   from  Eq.    (5.  lie),    i.e., 

With  the   premise   of 


*n  =    Ks   =    0, 

Eq.    (5. lie)   becomes 


when  written  in  terms  of  the  three  layers.  Application  of  the  funda- 
mental theorem  of  calculus  of  variations  gives 

o  Kta  i  Ypt 

s  °<3  =  -  ( —  +    * I  *  ^ 


*  [  J  Yj*  )  ?. 


o<3  °    *  U 
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c35_ 

iV 

0 

o  Y33 

•> 

6  (   ■=■ 

i( 

(5.49b) 


Here  we  come  to  the  well-known  but   interesting   paradox  of   thin  plate 
theory.      Note   that   the  variational   principle   does   not   yield  expres- 
sions  for      's  or     "S     ,    which   implies   that   they  are  negligible. 
However,    the   quantities      '  (f     and     "(f,  which  are   defined    in  terms   of 

S  and        S      ,      respectively,   appear   in  the  equations   of  motion, 

Eqs.    (5.19)   and    (5.20),    and   the  stress   boundary  conditions,   Eqs.    (5.37) 
and    (5.38).      In  addition,   the  general   stress-strain  equations   of 
Sect.   3.3.4,    in   particular   Eq.    (3.107),    show     that  the  condition 


'  33  a    33 

6       =    6       ■=    0 

is,  incompatible  with  the  basic   assumption  of   the   present  theory,   namely, 


// 


jr.-    -     t,.     -    0  . 


33  v  33 

A  complete  discussion  of  this  paradox  for  thin  plate  (and/or  shell) 
theory  is  given  by  Gol' denveizer  [25].  Here,  we  merely  observe  that 
above  irregularities  are  a  natural  development  from  the  variational 
principle  and  are  in  complete  agreement  with  the  fundamental  hypotheses 
of  thin  plate  theory.  Also,  it  is  now  apparent  that  of  the  two  sets 
of  equations  of  motion  (and  corresponding  stress  boundary  conditions) 
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the  higher  order  equations  are  preferable  because  the  quantities   'q 
and   Q   do  not  appear. 

Eqs.  (5.49)  give  the  stress  in  terms  of  the  strain  energy 
density.   For  the  deformation  theory  of  plasticity,  the  strain  energy 
density  was  established  in  Sect.  3.3.4.  for  a  three-dimensional  contin- 
uous medium.  Therefore,  the  results  of  Sect.  3.3.4  can  be  applied  to 
the  present  theory.   For  example,  one  form  for  the  strain  energy  den- 
sity is  given  by  Eq.  (3.98),  which  in  terms  of  subtensors  becomes 


'O 


(5.50) 


Partial  differentiation  of  Eq.  (5.50)  gives 

S9    (m       fit     6  vfi      **  LS    \ 


ivn     x 


(5.51) 


Because   of   the   symmetry  of      y.  .     and  Eqs.    (5.51),   the  stress-strain 
relations    for  the  middle    layer  are 
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l  3  X.     J        * 

3  =      Z    -=       i^  <5.52) 

Y 

o 

I22  =   2  *  Pj  *  Uak-zk  )  ?' 

where  Eqs.  (5..51)  have  been  substituted  into  Eqs.  (5.49).  Although  in 
a  different  form,  the  above  equations  are  the  same  as  Eqs.  (3.107). 
For  the  upper  and  lower  layers,  a  more  subtle  approach  must  be  taken 
because  of  the  inconsistencies  of  thin  plate  theory,  as  previously 
mentioned.   Eqs.  (5.49)  and  (5.51)  together  yield 

s    -  2  v  *    +3(3*''2*?)Yk  3    -      (5,53) 


The   condition 


5  33    =     0   3  (5.54) 

which   is   Love's   assumption,    is   applied  through  Eq.    (3.107)   and   gives 

i 

?  &  -  ?i 

i     3  ^  V         d  K  'W 

Y3      =       * /        .  (5.55) 

Next,  substitution  of  Eq.  (5.55)  into  Eq.  (5.53)  gives 
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S  l  =  2  P    Y  r  t  2  p  I  — \-  )    /     a  r   (5.56) 

r 

a 

which  reduces   to  the  expressions   for   classical   theory  as   given   in 
reference    [22], 

Similarly,    the  stress-strain  equation   for   the    lower    layer   is 

'i 

Y 

Also,  the  inverse  relations  can  be  written  immediately  by 
referring  to  the  equations  of  Sect.  3.3.U.   From  Eq.  (3.96),  and  with 


5      =    533    -    £  , 


the   inverse   relations   are 


*    'a-,  (5-58a> 
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S33  I   /    I  I     t   \   SV  I    /    I  I     \    -3 


'  -j(rrTj)^-j(rrl)^> 


V^       I    X   %<Q    .     L/l  _  1   X  \'c«     ■//> 


(5.58b) 


Naturally,  the  equations  can  be  written  in  different  forms;  for  example, 


S*  =    K^h    V. 


a1*-   i«n  r„  *  E^3yS3, 


-«3  -  c/3  Y3     -. 


-33  -  4833     -  -  3337      — 


(5.59) 


^  =   'V^s  V. 


rt       5 


where 


5. 


•^^  J»/  ^  ^     ^  ft\     ?  j/^-2  g'  \  ^  y* 
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5.2.   Buckling  Problem 

The  method  of  establishing  the  equations  associated  with 
buckling  has  been  developed  in  Sect.  3.4.  Therefore,  it  is  easy  to 
show  that  corresponding  to  Eqs.  (5.4)  the  displacement  functions 
Y.   associated  with  the  buckling  phenomenon  for  the  upper  and  lower 
layers  are,  respectively, 


'fi  =  -  (S3,i  *    I   $3,i)  , 


(5.62) 

Then,  the  infinitesimal  displacements  that  occur  during  buckling  for 
the  three  layers  are 

-f(bU3,i   *  |  h$a,i), 

hVL     =      hU-      t     9       bgi      ,  (5.63) 


iVL  -  bUL  -j('h$i*  S03}i  -\b$3,i) 
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The  variations  of  stress  resultants  have  been  established  and  are 
defined  by  Eqs .  (3.114)  -  (3.119);  however,  the  quantities   6el3  that 
appear  in  Eq.  (3.118)  have  to  be  established  for  the  current  theory. 

5.2.1.   Field  equations 

The   field  equations   are  determined   from  equilibrium  equations 
that   result   from  Eqs.    (5.26)   -   (5.29)   by  omitting  the   inertia  terms. 
The   first  equilibrium  equation  for  the   stable   position   is 

By  substituting  the  relations   between   the   stable   and   unstable   condi- 
tions,  e.g., 


into  Eq.    (5.64),   we   obtain 


Then,   by  subtracting   the  equilibrium  equation   for   the  unstable   condi- 
t  ion ,    i.e., 


N  ^+  N  r  +  N  '  )L  +  p  +  .p  +  p      =    0 


? 
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from   Eq.    (5.65),   we   obtain 

after  cancelling  the  6.  Eq.  (5.66)  is  the  first  of  four  field  equa- 
tions; the  remaining  equations  are  determined  by  the  identical  proce- 
dure,   and   can  be  written  by   inspection   from  Eqs.    (5.27)   -   (5.29) 


(5.67) 


r 


W*  n^  W<)  u„  HliV-Vh*']?., 

+  (sr?'  +/Y)P  +/#    J&391  +[  J  (M    -W    )  +  /MP  J  f3jf 

-  [  I  (j'-y<)  -  W-  j/i  +  >st  pi "/  -  o, 
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if  hi 

+  (i%-*")0+rs)+$(y->')+cJ-  o. 

(5.69) 
In  Eqs.  (5.68)  and  (5.69),  the  higher  order  terms  such  as 

have  been  neglected. 

To  complete  this  section,  the  quantities  P  and  C  must 
be  established.  First,  consider  the  term  'p  which  is  defined  by 
Eq.  (5.23);  for  the  stable  position,  the  corresponding  quantity  is 

/ 


zP 

Since 


,'o(  'q«3  'q-/3        <  >c/3  '{f=    1 


ft-*?*  if* 
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and 


AT'   6      =      ^       , 


we   obtain 


d)2r    =    d  ^    ,       at  u  -    j 


In  a  similar  manner,  we  obtain  from  Eq.  (5.2*0 


^     0      =     J      , 


Sj>>    «    S'e 


where 


be32-   's"(t>U3„  +  §  h§3l< 


<   '33 


h     r 


Summarizing,  we  define  the  following: 


S\ 


y  'c\ 

m         i  m 


>  *3         T,3 


=  M 


-43       -p.3    -°<3 
A,         a   JL 

a  43       y\3  »  43 


'p3 .  _<^ 


H 


P3   c3 

"p3    "c1 


)-H 


i   33  'J  i   33 

e     t  9    e 

-33  rf    -33 

a   33  n2  "  3J 

,  e   ,  6  e 


) 


\ 

u3 

'h 

( 9  y 

'i) 

)at{ 

f 

4 

"  3 

6 

"h 

J 

\      j 

Kz) 

(5.70) 
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which  are  used   in  conjunction  with  Eqs.    (3. 119).      In  the  above, 

(  i  33  i   I  d3  /    -  ,     h      ~  \  (    '  33 

he     =    6  a,    (  LL3yo<  +  -  y3)<  )  +    d  a. 


-33         c    -  <3  /    -  sJ   -  v        r  ,  33 

(5.71) 


h  e '  =  bsL    (  Us,*  +  9  y/j9at )  +  SJi,3  (/+  f3 ) 


<  "  33         c  ii  <43  f  -  h     -  \     ,      (  "  33 

be       -    dJb      (U.3,4  -    j  fj,j)   +    oa, 


5.2.2.   Boundary  conditions        < 

For  buckling,  the  boundary  conditions  are  obtained  from  the 
boundary  conditions  for  the  two  equilibrium  positions.  First,  consider 
the  stress  boundary  conditions.   For  the  top  surface,  the  equations  for 
the  3table  and  unstable  positions  (refer  to  Sect.  5.1.2)  are, 
respectively, 


2>*  -y  =  °,   y*-y  =  °>  *™ 

and 


y  -y  -  °,    y*  -y~°,  ^w 
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where  by  definition,  from  Eqs.  (3.49)  and  (3.110), 


/  ,  L 


i  *s  i  i 


P 


zr* 


zP 


i  r*s  (_  I  ~~>  (_  i 


-  t     =  t     +  be 


\      AT    d^^+'h 


or 


'  Irs     L 

-y  *  Kn , 

^ 

=y  +  hy  ■ 

Similarly, 

I       L 


'*  C 


'cL   =    'cc  +    &  c 


Therefore,  subtraction  of  Eqs.  (5.72b)  from  Eqs.  (5.72a)  yields 


1  n  c 


tni 


'  „  L 


'/*<> 


i-p;-  &t-r-  o,     i/;  -  s2c  =  0, 

which  are  the  stress  boundary  conditions  for  the  top  surface  during 
buckling. 

For  the  edge  surface,  the  stress  boundary  conditions  are 
obtained  from  Eqs.  (5.33)  -  (5.36),  (5.39)  and  (5.40).   For  example, 
in  the  stable  position,  Eq.  (5.33)  goes  over  to 


1U1 


o  o  o  o     A    \  / 


where 


&*  =  S^    SjJ 


0 


3 


V*«  V*  h'mJP     « 


■>     -c' 


Subtracting  the  equation  corresponding  to  Eq.  (5.73)  for  the  unstable 
position,  i.e.,  Eq.  (5.33),  from  Eq.  (5.73)  gives 

i  As  +  &  As  +   &A?=  m0  {  b'/w  "  +  h/Y)     +  S'/y?^  ) 

o  o  o  o     p    \  ' 

as  a  boundary  condition.  The  remaining  stress  boundary  conditions  are 
calculated  as  above;  in  summary,  the  stress  boundary  conditions  are  as 
follows  (for  convenience,  the   6  has  been  cancelled): 

on  the  top  surface, 

m 

on  the  bottom  surface, 

II  n  C  ll~,L  II      C  II      i 

P*  +   P   ■=  0         C     +  C '    =    0  ■  <5-75) 


1U2 


on  the  edge  surface, 


o^   +  t^       ^     ^0^6    \fflr+M+tf'j)  (5.76) 

li^'^J  +  cf    =  Mpl2{/*r-wp)  +  />n  rj  ,     (5.77) 

(5.78) 
J*      %  *       £ /*%*•*  \  -     ft     o(6      a     o(A        U     .,    JA     //46\1 


(5.79) 


Jp+£+U?  =  &{  (//*?+  N**+  i 


fi 


u 


3,4 


*[W-W)*tf]fM.+<!m%&+W)iL,ltl 


h    /    '   3       ii  3  v         7  J 

"7   (  ^  -^/   )  *  £       =    , 
Z    \.  0  0        /       o  o 


«f\li{**-'i<*h  **1  u. 


(5.80) 
3>*t 


'3,4    2lz\M  rtmr)-/mrt/yyir] 


(5.81) 
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In  the  above  equations,  the  definitions  of  Eqs.  (3.28)  -  (3.30), 
(3.114)  -  (3.117),  (3.119),  (3.144),  (3.145)  and  (5.70)  apply  to  the 
various  quantities  (i.e.,   'N   ,   'n   ,  etc.);  also,  all  higher  order 
terms  have  been  neglected. 

The  displacement  boundary  conditions  are  obtained  in  the  exact 
same  way  as  for  the  general  theory  (see  Sect.  3.4.2),  and  the  results 
are  as  follows: 

on  the  top  surface, 


V-     -       V ■      -       0   ;  (5.82a) 

on  the  bottom  surface, 


VL      -       VL     -  0-  (5.82b) 

and  on  the  edge  surface, 


i~> 


Vi  -    K  =  0, 


I/.  -       V  ■       =       0  (5.82c) 

L  L  1 

ii  „^ 

V;  -       V;      «      0. 


Again,    the    infinitesimal      6      quantity  has   been  cancelled   for  conven- 
ience. 


Wi 


5.2.3.      Variation   of   strain 

The   strain  variation   that   occurs   during  buckling   is   calculated 
by  the  method  of   Sect.   3.4.3.      From  Eqs.    (5.42),    the   strain  tensors   for 
the   stable   position  are 

-  '  •  *  JL 


for  the  three  layers.  By  substituting  the  equations  that  relate  the 
displacements  in  the  two  equilibrium  positions,  e.g., 

'it  -  'n  +  t%  , 

the  relations  between  the  strain  tensors  for  the  two  positions  are 
obtained  as  follows: 

> 
Yjo      ~        y^g      +        0  /J *      ~*       higher  order  terms; 


°<6        ~       ^*73  +       ^'o(A     "*       hi2her  order  terms; 
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o ^  -,      ~        0  /-  ~t"      0  /    ,  -      -t~      higher  order  terms ; 
«3j   =    '32,      ~*~     V  '   33       ~*~       higher  order  terms; 


//  // 


o(A  'g/A     ~*~      Of     j  a     ~t~      higher   order   terms ; 


where   the  strain  variations   are 


Ki  =  {(bv,>3^vh^vhjvM,hvM  «5„;, 

By  substituting  the  displacement  functions,  Eqs.  (5.7)  and  (5.63),  the 
above  become 


Ks  -  SA3  -  &JsA3  , 
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&rn  -  sr 


33 


33  ■) 


(5.84b) 


where 


*Kp-  kK?  -  ^\Kfi, 


fcfi'll  ty  *  VpU  +  |  (&L  -  fy<  -  Z  U3 


h«(S 


_  (5.87) 


(5.85) 


(5.86) 


(5.88) 


(5.89) 


(5.90) 
,fl<3      =     J    f/JW     t     f3     (p5j^     +     gi3    y)3itf    j    ?        (5<gi) 
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0^33     ~     f3    +    <p3   $3     ,  (5.92) 

+  {U3,J  ~   2   $3>4)(£3,p-  J   fj',/8)],         (5'93> 
//  - 

'    ^    ~    ~    (    Vsxtfi    '   J  p3,*{8)  >  (5.94) 

Note   that   the   strain  variations 

are    all    zero. 

5.2.4.      Variation  of   stress 

Consider   the  middle    layer   first.      In  Sect.    5.1.4,    all    the 
stress   components  were   found   to  be   finite;    furthermore,    the   stress- 
strain  equations   are    identically  the   same   as   for   the  general    theory, 
given   in  Sect.    3.3.4.     Therefore,   the   variation   of   stress    for  the  middle 
layer   is   also   identical  with  the  general   theory  as   developed    in  Sec.    3.4.4. 
Thus,    from  Eqs.    (3.162)   and. (3.163), 

Utf-     MCjU  Sru  (5.95, 
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where 

(5.96) 
or  in  subtensor  form,  Eqs.  (3.165)  are  applicable. 

Next,  consider  the  upper  layer;  two  methods  may  be  used  to 
establish  the  stress  variation.  One  approach  is  to  use  the  stress- 
strain  equation,  Eq.  (5.59),  of  Sect.  5.1.4,  i.e., 

and  operate  upon  it.  The  second  method  is  to  use  the  general  equations 
of  Sect.  3.4.4  and  reduce  them  for  this  special  case.  The  second 
method  is  more  general  and,  therefore,  will  be  employed;  the  equations 
that  will  be  used  are  Eqs.  (3.165). 

Recall  that  the  stress-strain  equation  for  the  upper  layer  was 

33 
developed  from  the  condition  that   7S   =0,  which  implies  that 


b  A,        -     0.  (5.97) 

Substitution  of   Eq.    (5.97)    into  Eq.    (3.165)   gives    (the   primes   are 
omitted   for   convenience) 


Mn"irti+2Mnair„  +  Ar's/'sl  -  o. 

3333 
From  Eqs.  (3.166)  and  (3.108),  it  is  seen  that  M      is  a  scalar 

quantity,  namely, 
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33 
where     S   '    =   0     has   been  used.     Therefore 


1333 

M3333   a/  n       MlsT3   M't3  .  <5-w 


,    .    I  O  J  J 

M  <  n  M 


01 33  M3323      61  n  jp^     h/t 

Substitution  o£   the   above    into  Bqs.    (3.165)   gives 


I  M3333  J    0/  ^ 

{  M3233  )°'13 

\  .  ^j  3333  J    v  I   t 


(5.99) 


After  some  lengthy  but  straightforward  calculations  the  above  can  be 


written   as 


Ml 

0/  n     +     Z/K    '  d/-/3  , 

(5.100) 


U-/rlhrn  +  2K*n*rnt 


W-  /r*irn  *  m™  irn 
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where 


R      J  =  A         ■>  B      J  3  (5'I01) 

A     =J(99  +99  )^1A^~1)9  9J> 

K  (5.102) 

M  _  A 

,      U       K 

(5.104) 
The  above  are  the  general  stress  variation  equations  under  the  conditions 


33  f       33 

6         =  &A,        «  0 


Note   the   symmetry  relations 

«<.ii  ^iotX'i  *L'iY  Y,' all 

R         J    =     R         '      -     A7      ^       =     /?    '  .  (5.105) 


For  the  more  restrictive  theory  under  consideration,  the 
conditions 
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i  o<3  ' 

-5  =  y<3  =  o, 


must  be  applied  to  the  above  equations.  Thus,  for  the  upper  layer,  the 
stress  variations  are 

Similarly,  for  the  lower  layer, 

When  reduced  to  classical  theory,  the  above  agrees  with  the  equations 
given  by  Teodosiadis,  Langhaar,  and  Smith  [22]. 


5.3.   Modified  Buckling  Problem 

As  discussed  in  Sect.  3.5,  the  modified  buckling  problem  is 
formulated  by  neglecting  the  effects  of  the  displacements  prior  to 
buckling.  Therefore,  the  equations  of  Sect.  5.2  are  applicable  when 
the  displacement  functions  u_   and   f,   are  omitted;  note  that  this 
is  the  same  as  assuming  that  the  plate  is  flat  prior  to  buckling. 
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5.3.1.   Field  equations 

In  view  of  the  above,  the  field  equations  as  obtained  from 
Sect.  5.2.1  are  as  follows: 


(5.109) 


'P 

-ie^->"h'c^c^n^f}L 

'    (5.110) 
-\  U.  f '   rf    "  <<fi\       '    °<6      "    c<A~lt  '-,3      -3     i/3 


+  ~t('ps-V)  +  cs  •  o. 

33 

In   addition,   the      6e   '       quantities,   which  are   employed   in   the  defin- 

.     .  3 

LtLons   of      6nP   ,      are   obtained   from  Eqs.    (5.71)   and  become 

S'C       =     Skf     +    S        (  £>L/3)0(    +    |     b  <f> 3)0(    )  j  (5.112a) 
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(5.112b) 

5.3.2.   Boundary  conditions 

By  neglecting  the  u3   and  73   terms  from  the  equations  of 
Sect.  5.2.2,  the  stress  boundary  conditions  on  the  edge  surface  become 

#  +  ?+#  -.*A  CmV+tf'+W), 

(5.115) 


/  "/  "  I  ($  +P)  VJ,  fa*'-  W- *(M*+tf»)], 


h  '#-#)  ./>  .*  {[(iftw+  i+)  *  r«]i 


(5.118) 
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The  equations  for  the  stress  boundary  conditions  on  the  top  and  bottom 

surfaces  are  the  same  as  Eqs .  (5.74)  and  (5.75)  of  Sect.  5.2.2.   How- 

33 
ever,  recall  that  the   6e  ~   quantities  have  been  redefined  per 

Eqs.  (5.112). 

The  displacement  boundary  conditions  are  given  by  Eqs.  (5.82). 


5.3.3.   Variation  of  strain 

By  reducing  the  equations  of  Sect.  5.2.3,  the  strain  variations 
are  as  follows: 

(5.119) 
(5.120) 

^i  =  I  (&&  +  S&3,<  +  9  Sips,*)  ■■>        <5-l21> 

sr„  =  sf3  j  (5.122) 

*  *  ifr,^  )  -  2B3(hU3}^  I  6<p3}y  )]. 

(5.123) 
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5.3.4  Variation  of  stress 

No  change  occurs;  the  stress  variations  are  given  by  Eqs.  (5.95), 
(5.106)  and  (5.107). 


6.   NUMERICAL  EXAMPLE 


Let  us  now  consider  a  simple  but  practical  problem,  namely, 
to  determine  the  critical  buckling  load  of  a  simply  supported  rectan- 
gular sandwich  plate  subjected  to  a  uniaxial  compressive  load  of  suf- 
ficient magnitude  to  induce  stresses  beyond  the  proportional  limit. 
The  plate  is  shown  in  Fig.  3. 
.3 


Fig.  3.   Rectangular  sandwich  plate 
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Because  of  the  geometry,  the  curvilinear  coordinate  system  becomes  a 
rectangular  Cartesian  coordinate  system.   Therefore,  all  contravar- 
iant  tensors  become  Cartesian  tensors,  covariant  derivatives  become 
partial  derivatives,  metric  tensors  become  Kronecker  deltas,  etc. 

To  obtain  a  solution,  the  equations  of  Sect.  5  will  be  employed; 
in  particular,  the  modified  buckling  equations  of  Sect.  5.3.   In  addi- 
tion to  the  assumptions  of  Sect.  5,  the  analysis  will  initially  be 
based  on  the  following  two  assumptions: 

1.  The  face  layers  are  very  thin,  so  thin  that  they  have 
no  resistance  to  bending. 

2.  The  middle  layer  is  a  so-called  soft  core;  that  is,  the 
only  significant  core  stresses  are  the  transverse  shear 
stresses . 

The  first  assumption  requires  that 
'h 

„       Xp    QZ      J03       =  0.  (6-2) 

I 
Recalling  Eq.  (3.3),  i.e., 

Eq.  (6.1)  can  be  written  as 


[   '%p[*3-{(h+l)]J§3=   0. 
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From   the   definitions   of   Eqs.    (3.28),   the  above  becomes 

My    =      i    (h  +    %)   A/<fi    .  (6.3) 

Similarly,    through  Eqs.    (3.2)   and    (3.28),    Eq.    (6.2)   can  be  written   as 


My    =    ~    {    (h  +  "h)    Ny.  (6.4) 

m 

Obviously,   equivalent  relations   hold   for   the  variational   stress 
resultants;    then, 


S/my  =   {  (  h+'h)  Shy   ,  (6-5> 

t>""y      =    "    {     (*>+'*>)    i%y    .  (6.6) 

The   second   assumption   requires   that 


^33     =      Sy     -      0,  (6-7) 

which   implies   that 


dA^33     =       O^L^a    ~    0  .  (6.8) 


Therefore,  all  the  stress  resultants  that  depend  upon  the  above 
stresses  are  also  zero,  e.g., 


(6.9) 
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Also,  in  Sect.  5.1.4,  it  was  found  that  under  the  premise  of 


Y-3  -   0 


the   variational   principle   yielded  the   conclusion   that 


3C3     ~      &'t3      '       ^^L3      =       0. 


Similarly,   the   assumption  of   Eq.    (6.7)    implies   not   only  Eq.    (6.8) 
but   also 

y33  -  r^  =  irn  =  hF^  =  o.  (6.io) 

When   substituted   into  Eqs.    (5.46)   and    (5.92),   the   above   yields 


^3     =      <p3     =     0. 


(6.11) 


When      \|[,    =   0      is   substituted    into  Eqs.    (5.4),   we   obtain 


=      -   LL  ,      ,      ;  (6.12) 


jfV   =    ^    =   -  ^w   i 


and,  thus,  the  equations  relating  the  displacement  functions,  corre- 
sponding to  Eqs.  (5.6)  and  (5.7),  are  as  follows: 


(6.13) 


"l 
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-      /  f  '  "  I    /  'i       "i    \    '  1  ~  l   /  '        "      \ 


The  equations  for  the  displacements  of  the  three  layers  become 


3   5 


(6.15) 


Similarly, 


3  ) 


4«4"lA*I  ^JW  >        ^  ■  ^3  ,    «.I6) 


and 
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6.1.      Field   Equations 

In  Sect.    5.3.1,    four   field   equations  were  established;   however, 
of    this   set,    Eq.    (5.111)   was   a  result   of   the  arbitrary   quantity     6?3 
in   the   variational   equation.      Because     "*3      is   zero  for  this   problem, 
Eq.    (5.111)   must  be   omitted.      Then,    with  Eqs.    (6.5),    (6.6),    (6.9)    and 
(6.11)   substituted    into  Eqs.    (5.108)    -    (5.110),    the    field   equations 
become 

(6.20) 
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Next,  by  assuming  that  the  variation  of  body  forces  are  negligible  and 
by  applying  stress  continuity  at  the  interfaces,  the  definitions  of 
Eqs.  (3.119)  and  (5.70)  along  with  Eqs.  (5.112)  and  (6.11)  yield 

p.  ,  p.  .  P.   =  A  _  A    t 

With  the   above   and  by  using  the   first   two  field  equations   to  simplify 
the   third   equation,    the   field   equations   become 

blfrp     *     ^0  7/9      +    ZP<X    "fat        ~        01  <6'22) 


+  fc-fc-0;      . 


where 


-A 


h    [  Z 
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6.2.      Boundary  Conditions 

The   basic  equations   for  the   boundary  conditions  were  established 
in  Sect.    5.3.2.     As  mentioned   in  the   previous   section,    the   fact   that 
6ijr,      does   not  exist   for  this   case   requires   the   omission  of  certain   equa- 
tions,  namely,   Eqs.    (5.116)   and    (5.118).      Substitution  of   Eqs.    (6.5), 
(6.6),    (6.9),    (6.11)   and    (6.21)    into  the   remaining   equations    yields    the 
stress   boundary  conditions   as   follows: 

on   the   top  surface, 


'p  ■  -  F '•  =  o       'c  ■  -  'c  ■  -  0  ■ 


on  the  bottom  surface-, 


on   the  edge   surface, 


(6.26) 


(6.27) 


A  *  £<  +  *<  "  ^  (^  +  Kfi  )  ;  (6'28) 

jd-^w-A^h/w  m^C^f^f)z '  (6,29) 

A^-I^AJ-M^v*^3,  (6*30) 

+  ?(h<n<f'h*<fi)*  +h'fr  Vffil-  (6'31) 
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The  displacement  boundary  conditions  are  as  follows: 
on  the  top  surface, 


]/.    -  ]/.     =  0  (6.32) 

on  the  bottom  surface, 


y.    -  /•  =  0   .  (6.33) 

on  the  edge  surface, 


V;   -    V;    =  0 


y.   _  V )    =  0  ,  (6'34) 

Y  L  V  L 


6.3.   Variation  of  Strain 

With  Eqs.  (6.10)  and  (6.11),  the  strain  variation  equations 
of  Sect.  5.3.3  reduce  to  the  following: 
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^3   =  Z    i^f*     *      t>V3,«)   )  «.36) 

+  zWa,*p)-*B,ida,«].  .-37, 


6.4.   Variation  of  Stress 


With  reference  to  Sect.  5.3.4  and  with  Eq.  (6.8),  the  stress 
variation  relations  are  given  by  Sqs.  (5.95),  (5.106)  and  (5.107); 
by  substituting  into  these  equations  the  strain  variations  of  Sect.  6.3, 


we   obtain 


SJL^      =       MJU3     (Sfr    f     hUSti)  ,  (6.38) 

*  2U3,n)  -   203  SU3,ri  ]. 
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Because  of  the  symmetry  relations  of  Eq.  (5.105),  we  have 

R<pn  hU^h    =    ^/srs   hU^Y    ,  etc* 

Therefore,  the  stress  variations  can  be  written  more  compactly  as 
foil ows : 


&<**  '  Kant^ra  + 1  (ift,i  +MJMi)-B3  lU,,n], 


(6.39) 


S^3     =    Mjm     (8(pY     +     6l73,y)}  (6.*0) 


(6.41) 


6.5.      Solution 

As  previously  mentioned,  the  purpose  of  the  modified  buckling 
equations  is  to  avoid  solving  the  prebuckling  problem.   However,  for 
a  formal  solution,  the  initial  state  of  stress  in  the  sandwich  plate 
should  be  determined  from  the  prebuckling  stress  boundary  conditions; 
in  particular,  Eq.  (5.33)  which  reduces  to 
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for  this   case.      If   the  uniaxial   compressive    load,   say     -  R   ,    is   taken 
to  act   in   the      92     direction;    then,   Eq.    (6.42)   yields 

N2Z    +    Nzz      =      -    Rz     7  (6.43) 

k,    +    V//       »       Nn    +    A/,2     -    0  .         '  (6.44) 

During  buckling,  the  upper  and  lower  surfaces  are  free  of  tractions, 
i.e. , 


this  is  a  boundary  condition.  Therefore,  Eqs.  (6.26)  and  (6.2?)  yield 


/ 


A  -  A  =  o 


z'l      ~    fri      =  U.  (6.45) 

The  remaining  surfaces  for  which  boundary  conditions  must  be  specified 
are  the  edges.  Because  the  edges  are  simply  supported,  it  is  known 
that 

which  from  Eqs.  (6.34)  give 


K  =  V.     =  V,  -  0.  (6.46) 


As  for  the  other  displacement  components  on  the  edges,  we  only  know 
that  a  motion  in  the   P   direction  must  occur  at  the  edges 
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9^=0,   L.   All  other  displacements  may  or  may  not  be  restrained; 
fortunately,  we  can  leave  this  decision  unresolved  for  now. 

From  Eqs.  (6.43)  -  (6.45),  the  field  equations  become 


—  ^ 

-    R2     U3722     +    H     5         ^        =       0.  (6.49) 

Recalling   the   definitions   of   the  stress   resultant  variations,   e.g., 

a, 


i>Kp  "\b   h^^    Jh   , 


the  stress  variations,  Eqs.  (6.39),  can  be  applied  to  obtain 

a. 


Integration  yields 

Similarly, 

h     7  % 


K/s  *  h  Ryu  (  uy,i  "zf^fi  u3>  n  )> 


(6.51) 
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fjt  =   f>  Max*  (Pr  +   U»t)-  (6'52) 

After  substitution  of  Eqs.  (6.50)  -  (6.52)  into  Eqs.  (6.47)  -  (6.49), 
the  equations  become 


+  h  ,7.  .        -  'h  H 

■>Y6f  J 

(6.53) 

h     r.  %   n 


ri     T  h    ,7 


-2 Kan  (h  +  ^.^    =    ^, 


<\    -       — 


-  R2  Us,zz  .■*■  h  h  M<M  (tpr,<  *  0S,U  )  =   0. 

(6.55) 
To  obtain  a  solution  to  the  above  equations,  concepts  employed  in  the 

work  of  Chang,  Ebcioglu,  and  Baltes  [26]  will  be  applied. 

It  may  be  recalled  that  in  Sect.  3.2  it  was  mentioned  that  the 
basic  equations  could  be  formulated  in  terms  of  two  sets  of  displace- 
ment functions.   Examination  of  Eqs.  (6.16)  shows  that  for  this  problem 
the  field  equations  can  be  simplified  by  converting  to  the  alternate 
set  of  displacement  functions.  Thus,  substitution  of  Eqs.  (6.16)  into 
Eqs.  (6.53)  -  (6.55)  gives 
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/,     /_  /  ,  //,     //_  // 


"ZMlrtiljCl/y-Uj+Hht]*     0,  (6.57) 

-  R2  'VM2*h~h  M<3„[j(i/f„-  (/,,<)*%  U„u]'  0. 

(6.58) 
Eq.    (6.56)    is    identically  satisfied    if 

/.        I  II        u 

"   °Jp*i    =     ^    h   ^«/3*6   ;  (6-59) 

Uf    =    -    LU     Uf   j  (6.60) 

where     uu      is   a   constant.     Therefore,    substitution  of   Eqs .    (6.59)   and 
(6.60)    into  Eqs.    (6.56)   -    (5.58)   yields 

0-  0, 

- Rt  i/,tU  +  h  h MM  ( '-f  'Uu  f  h  U3l u)-0. 

ft  t(L      Sl^ 


It  has  been  established  that  the  only  significant  prebuckling  stresses 
are 
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/  //  —  — 

Then,    Eqs.    (5.101)    -    (5.103)   show  that 


—i  ' 

/v//2  ^ 1121  /2U  **2IH  U 


/—i  '  ' 

K J  222   '      '^2122  ^2212  ^2221     ~    &\ 

therefore,  when  expanded,  Eqs.  (6.61)  and  (6.62)  become 


/       r   ' 


-Mm  (-J-  U.  +  mJ-Mmf-J-  4^4,*M 

(6.63) 


H  {Rm  *  R<2/2  )  U/,tZ   +  R/2J2  Uzrt  *  ^2222   U2>2l] 

(6.64) 

£0+u))[Mmi/„,  +Mm  (tihz  +  k,i)  +  M2323'{J2,z] 

+  h  h  (Mt3l3U3,n    *  ZM/323'(J3)IZ)  (6.65)- 

r(fi  h  M2323  -  Rz  )  y3       =  0, 
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where  (the  primes  and  bars  are  omitted  for  convenience) 

Ml323  =  rp(T7  ~  "tM^  , 


""a"       Y.       .?.<;2l  SY.        Y,  K*"'    ' 


(6.67) 


(6.68) 


o  "0 


2 


JA        r      Y0  (6.69) 

JK       7     t  (6.70) 

Rail  '-  T  (6'71) 


^-4*HK^M^**m3  . 


(6.72) 


From  Eqs.    (3.91)    and    (5.104) 


4  = 

J    °22   > 

/ 

^  (  1% 

v. 

(6.73) 


0=T 


(6.74) 


S,  w   -,      ,  '4      A 


^^.,Aj(^+JA+^J 
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Although  a  simplified  large  strain  theory  is  being  employed, 
it  is  usually  assumed  in  the  literature  that  the  stress  and  strain 
tensors  can  be  interpreted  in  the  classical  sense.  Then,  the  material 
constants  of  classical  theory  are  applicable,  e.g., 

r  _         E 

6  ~  ITJ^J  > 


k.  = 


(6.76) 


3  0-Z/l)    ' 

where     G      is   the  shear  modulus,   E      is   Young's  modulus,      (i     is   Poisson's 
ratio,   and      k        is   the  bulk  modulus.     As   shown  by  Teodosiadis,    Langhaar, 
and  Smith   [22],    the   following   relations   hold: 

r.  '  6*  =  3-0-2^)$.  >  <6-77) 

iS>    =    Q =     £* .  (6.78) 


where     G        and     G       are   the   secant   and   tangent  shear  moduli,    respec- 
tively,  and     E       and  E        are   the  secant  and   tangent  moduli,    respec- 
s  t 

tively,    as   obtained   from  a   simple   compression  stress-strain  curve. 

In    1950,   Hoff    [27]   found   the   solution   to  the  elastic   buckling 
problem  of   a  rectangular  sandwich  plate  by  assuming  displacement  func- 
tions   as    follows : 


17k 


U,  -  J,    cos  iFj —   s/At  — —    • 


',  i        V  air '6,  nnTT  69 

Uy  ~    h    SIN   -£-7—      COS    — •  (6.79) 


™  , 


',  -       v  air  9.  /n  n  9Z 

(J3  =    J3     S/N  f-y-'     SiH   —J—     . 


These   same   relations   were   used   by  Chang,   Ebcioglu,    and  Baltes    [26] 

to  solve   the   plastic  buckling  problem.      The   above  displacements   satisfy 

the   simply  supported   conditions,    i.e., 


U3  =  0}        AT     BrA[)      ANO     S2- 


'  0    I 


and   also   impose   the    limitations: 


'Uz-0,     tr  6,41). 


Substitution  of  Eqs.  (6.79)  into  Eqs.  (6.63)  -  (6.65)  gives 


Actually,  q  was  taken  to  be  one.   Here,  we  allow  for  higher 
buckling  modes. 
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V  \2      'l    'n         /  MTT  \Z       i 


^KAH-J^KAtJ+'-it^X 


i.    ,  / 


-/- 


H(%)MmY3)co5$¥>SM 


m 


rr9? 


L 


1 


Y,  +  t 


h 


L 


V 

(6.80) 


{Ww,+  kim)(J£x!r)X 


L>  A  Lt 


*  o>^m  *^m{%p  ^  a„u 


h 


1T0Z 


\l(i^){$t)Mm  Y,  +  t  (/tu,)(f)  M23Z3  Y, 


(6.81) 


£  _ 


+ 


+ 


P  Hf)  fad  hta-/aff]rs}*»*P°»al& 

t)Ml*3Ys]co5lf>caS<fe>~0. 


(6.82) 
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If     M._„»      is   zero,    the  equations   can  be  written  as 


d,z  d>/3      \    f     X     \ 

d21        dz3      \\      Y2     I    =    0 

^32  ^33    I    \       '3 


where 


J  'r,         2       'D        ("lL*\Z       G  (l+uj)  /   L*\ 

dn  "  A////  f  +  Kan  \  Lu)         'hf,       I   ir/ 

a 


2 

r  I     ? 


d/2    -    d2,    ~      (  /?„2Z    i-    Kml  )   V1 


M  L 


K 


i      _     o  h  &  L« 

cL/3   -         7TD        , 
If  h  1 


-a 


^22     *    A/2/2  f    +   "2222  [      /      /T         ^  £  \     fT  *       1 

d33  ~  G  ^>  h    ^-  +  (G  h  h   -  f?2)(  —  )    ; 

and   a  unique   solution  exists    if   the  determinant   of   the     d.  .     coef fi- 
ll 

cients    is   zero.      For     M,__«      to  be   zero,   Eq.    (6.67)   requires   that 
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hS0         60 


or       G ,<  -    6 


if.      I   '  *' 

which  implies  a  Hookean  material,  which  was  assumed  in  reference  [26], 
Then, 

S  -  6S  =  St   ; 

therefore, 

The  middle  layer  is  elastic. 

Because  R   and  the  parameters   'r      depend  upon  the  state 
of  stress  prior  to  buckling,  an  explicit  expression  for  the  critical 
buckling  load  cannot  be  obtained.  However,  parametric  curves  can  be 
established;  following  the  procedure  of  reference  [26],  the  following 
parameters  are  introduced: 

p  -    IEAL)         ^  E  h  E 


L,    I  {l-'/SXt+mt)    '    /m'    %  E      ' 

£   /  1  -f    -   -A 
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£ 

»           F° 

6    - 

a-y)  'h    > 

1  ~rie  > 

Z." 

Run 

'H 

-7                 Kj/22 

Li2  '        '#         7 

V 

R2221 
'ti 

-7                1     f^/2/2 

t33                 <H            , 

(6.83b) 


where   (R-)    is  the  critical  stress  resultant.  With  the  above  para- 
2  cr 

meters  substituted  into  the  equations  for  the  d.  .  coefficients,  new 

L] 

coefficients  are  established  as  follows: 


7         2..  2  €  p 


7 


2<i.    j  -  Z. 


fa  -ft*  -  -f  (z«+  ~i3)> 


ir  h  &  j  e 


'13 


f6 

h3        r  A  /&    ; 

A  =  ^  * 
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Then,    after   some   rather   lengthy  algebraic  manipulations   the  equation 


f?/i  fiz  fit  3 

Pzi  A?2  Hz  3 

1^31  P32  @33 


0 


yields 


Kp 


-  2ZIZ  ?3i  ^  p  -  fz*  f  f  +Z„  Z33fp  )rA 


+  Bfefk  (Zzi  +  **ff)]-  ZfY€f*  ( 


2, 


u 


+ 


?22-?33-2?iz)}x{rp[  i 


Z-22    ^33 


2       2 


f  1  Z„  Zzz  J  p   ~    2in  Zxjf  -  iZlz  j 


z    z 


+  I„Z 


■f 


When  the  equations 


/*/*  [i  (f+jY)z«  rp< + *Y*fi  $ 


(6.81) 


*J« 


I      '  H     II 


h  S22+    In  527 


(6.85) 
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D   "    (flh)       ~ 7 ,  (6.86) 

h  +   n 


f 


j 
"c         ; 

^  21 


(6.87) 


are    introduced    into  the   definition  of     K    ,    i.e., 

P 


Kn  = 


m 


Cl? 


'  o 
a  rigidity  parameter   i3   obtained,   namely, 


2 ir  D  (%+%)  0  ,     ,  2]  (Sul 


2   (/-/t    J   — .       (6.88) 


'h('h  +  l)(^)lf  K    ^'     *P 

7 
Eqs.  (6.84)  and  (6.88)  constitute  the  solution  to  the  problem.  For 

a  given  sandwich  plate  under  a  given  state  of  stress  in  the  facings, 
the  buckling  coefficient  K   can  be  calculated  from  Eq.  (6.84)  for 
various  buckling  modes,  and  the  rigidity  parameter  can  then  be  deter- 
mined from  Eq.  (6.88).   Comparison  of  the  calculated  rigidity  para- 
meter with  the  actual  plate  shows  if  buckling  will  occur. 

As  a  check  on  the  validity  of  the  analysis,  the  equations  will 
now  be  compared  with  the  work  of  Chang,  Ebcioglu,  and  Baltes  [26]   in 
which  the  following  was  assumed, 

1.  The  face  layers  are  the  same  material. 

2.  The  state  of  stress  in  the  face  layers  is  the  same. 

3.  The  middle  layer  is  orthotropic. 
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4.  The  plate  buckles  in  the  simplest  mode. 

5.  The  face  layer  material  is  incompressible. 

6.  A  correction  factor  can  be  applied  to  account  for  com- 
pressibility as  given  by  Seide  and  Stowell  [28], 

! 
From  assumption  #1  and  Eq.  (6.83), 

%     ■ 

<W    =  ~TT      •  (6.89) 

h 

From  assumption  #2  and  Eq.  (6.87), 

=  /.  0  .  (6.90) 

From  assumptions  #1  and  #2  together, 


and,  therefore,  from  Eq.  (6.59), 

/ 

h 

LU     ~     T-  .  (6.91) 

n 

From  assumption  #3,  reference  [26]  has  defined  the  parameter 


? 


6 


X 


s 


(6.92) 


which   is   the  ratio  of   the  core   shear  moduli   in   the      9.    and    P       direc- 
tions,   respectively.      For  comparison,    it    is   necessary  to  take 


3 


~      1.0.  (6.93) 
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Assumption  #4  requires  that 


f 


-1.0  <6-9^ 


see  Eqs.    (6.79).     Assumption  #5   implies    that      'u.  =    1/2.     Therefore, 
the  material   constants,    from  Eqs.    (6.76)    -    (6.78)   become 

/C,  =  °°,         &s=    j£5    ,  G{    «   J    £{     .  (6.95) 

Then,    from  Eqs.    (6.69)   -    (6.7U)   and    (6.83), 


W/         £33  f-Oj  z/Z  2_    -> 

(6.96) 

From  assumption  #6,    the   corrected   critical   stress    is   given  by  the 
equation 


'/< 


z 


where  K   is  the  elastic  buckling  coefficient  given  in  reference  [29] 


as 
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For   the  special   case   of      g'  =   1.0,      the   above  reduces   to 

Kt  -  (i-t-p2)    l(H{5z )r  +  f$  j  (6.98) 

which   is    independent  of  Poisson's   ratio.        Therefore,   Eq.    (6.97)   becomes 

3_ 

When  Eqs.  (6.89)  -  (6.91),  (6.9U)  and  (6.96)  are  substituted  into 
Eqs.  (6. 84)  and  (6.88)  the  result  is 

Kp  =  ^{(/^)[i(£2^¥Z„/32-2/6^/3v)r  §■ 

+  il2Zf>1-  2/3* + {3* ) r  $  +  2f,2(Z2J*{32)] 
-i-  f>Z  I    J  (l*ji')r$  f    2/32]j       ,        (6.100) 


2ir2D       _  3_    feJ^  (6.l01, 

Ll(U%)    ~  Z       KP        ' 

which   is    identical  with  the  equations   of   reference    [26]   for     g  =    1.0. 
Therefore,   the   solutions   agree. 

As  mentioned  by  Nadai    [30],    Poisson's   ratio   is   variable   in  the 
transition  region  of  elastic   to  plastic  stress,   which  suggests   a  pos- 
sible  refinement   for  the   correction   formula  Eq.    (6.99),   namely, 
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(U  =11^^ 


f  I  ~c*lp.-  z 


where      jj,        is   the  variable  Poisson'3   ratio.      Nadai   [30]   gives   an  equa- 
tion for  the   approximation  of      p,       as 


v 


Y 

where      y,       is   Poisson's   ratio   in  the   plastic   range,      y        is   the   elastic 
strain  at   the   yield   point,   and     y  is   the   total  uniaxial   strain. 

For  uniaxial   stress,   Eq.    (3.96)   can  be   solved   to  show  that 

Y   =-  $ 
7       £    7  ' 

/  =  —  c        /  >     y 

where  S    is  the  yield  stress  in  compression.  Then,  with  the  usual 
assumption  that   u.  =  1/2,  Eq.  (6.103)  becomes 

(6.104) 

It  should  be  realized  that  the  expressions  for  K   and  Z,,, 

P        11 

Z   ,  etc.,  are  simplified  by  the  assumption  of  an  incompressible  face 
material;  however,  this  assumption  can  readily  be  waived  if  one  has 
access  to  an  electronic  computer.  Therefore,  it  should  be  of  interest 
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to  compare  the  results  obtained  with  and  without  assumptions  #5  and 
#6.   For  a  compressible  face  material,  Eqs.  (6.84)  and  (6.88)  become 

KP  =  {('*?)[*  {Za  233  *  <tZ„  Zuj3Z-  ZZa  luf 

-  1il(?  * 1„  z»f)r  +  2(f&  (in  i-  Z„  (?  )] 

-  sfe  (i„  1 4  -  Z„  -  2ZU  )}<{r[i  (Za 1. 


33 

z     z 


+  ¥£M  I2Zp  -  ZZIZ  Z33/6  -fZ/2/6   +zlt  z33p  )r 


2rr2  D      =  Zd-'ic2)  — ^ 
Lx(h+h)  kp 


-     Z(i-)L2)    ~~     ;  (6.106) 


when  Eqs.  (6.89)  -  (6.91)  and  (6.94)  are  employed.   In  these  latter 
equations,  the  general  expressions  for  Z   ,  Z   ,  etc.  are  to  be  used, 
i.e.,  Eqs.  (6.69)  -  (6.74)  and  (6.83). 

To  obtain  actual  numerical  results,  a  computer  program 
(Fortran  IV)  was  written.  The  program  along  with  a  description  of  its 
usage  is  given  in  Appendix  A. 
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6.6.   Results  and  Conclusions 

For  analysis,  a  square  sandwich  plate  (actually,  3  =  1.0) 
with  17-7  PH  stainless  steel  facings  was  chosen.  The  properties  of 
17-7  PH  stainless  steel  are  given  in  Appendix  B.   Under  the  assumption 
of  both  compressible  and  incompressible  material,  solutions  were 
obtained  for  a  wide  range  of  stresses,  extending  well  beyond  the  0.02% 
yield  point,  by  using  the  computer  program  given  in  Appendix  A.  The 
results  are  shown  in  Figs,  k   and  5. 

Fig.  U   clearly  shows  that  if  the  simpler,  incompressible  the- 
ory is  used  the  results  must  be  corrected.  Although  conservative,  the 
correction  given  by  Seide  and  Stowell  [28]  deviates  considerably  from 
the  correct  (compressible  theory)  solution  at  the  higher  stresses. 
Much  better  agreement  is  obtained  by  using  the  variable  Poisson's  ratio 
in  the  correction  equation;  unfortunately,  however,  this  procedure  pro- 
duces nonconservative  results  for  an  r  >  .5,  where  r   is  the  core 
shear  stiffness  parameter. 

The  solution  for  compressible  theory  is  shown  in  Fig.  5.  With 
r  as  a  parameter,  several  curves  are  shown. 
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APPENDICES 


A.   COMPUTER  PROGRAM 

The  solutions,  as  given  by  Eqs.  (6.99)  -  (6.102)  and  (6.10*0  - 
(6.106)  are  evaluated  by  means  of  the  accompanying  Fortran  IV  computer 
program.   The  input  data  required  are  the  material  properties  of  the 
face  layers;  in  particular,  each  data  card  must  have  the  following 
information:   compressive  stress;  Young's  modulus;  tangent  modulus; 
secant  modulus;  and  Poisson's  ratio.  The  data  are  to  be  punched  under 
the  format  (4E10.0,  F10.0).   Naturally,  as  many  data  cards  can  be 
employed  as  desired.  With  this  in  mind,  the  first  card  of  the  data 
deck  must  specify  the  number  of  following  cards;  and,  in  addition, 
the  0.02  per  cent  yield  stress.  The  format  for  punching  the  first 
card  is  (14,  E10.0). 

Besides  the  input  data,  the  equations  are  dependent  upon  the 
parameters   6  and  r,  which  are  controlled  by  program  deck  cards. 
If  an  aspect  ratio   0   other  than  one  is  desired,  say  four,  then  the 
card 

B  =  1.0 
must  be  replaced  by  the  card 

B  =  4.0 
The  core  shear  stiffness  parameter  r   is  controlled  by  the  two  cards 

70  R  =  0. 

50  IF  (R-l.l)  40,40,20 

The   output   is   completely  self-explanatory  except   that  the    letter 

P      is   used   to  represent  Poisson's   ratio.      The   program   is   as    follows: 
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PLASTIC    B  JC    LING 
ITE     (     ,     10) 

10  FORMAT  (72H  PLASTIC  BUCKLING  if  IICH  PLATE  \ri  I  TH  COMPRESSIVE  LG 
IAD  IN  Y  i  .  /  661  ER  OF  HALI  AVESt  IN  X  DIRECTION  =  i.O 
1,        IKY,   II      ..I  i     N    =     1.0) 

-     i  . 
WRI TL     (b ,      :      ) 
»0    FORMAT     (15H    ASPECT     tATIC     =     ,    (-3-1) 

15,    65)    J i 
65    FORMAT     (14,     E10.0) 

ui  i  -  ,    ,  J 

READ  (5,30)  S,   ,  [J  I  ,  DS,  P 
30  FORMAT  KE1  .Oi  F10.  ) 
(b,i,l 

11  FORMAT  (5X,  >8H  COMPRESSIVE  olRESS  YOUNG'S  MODULUS  SECANT  MODULU 
IS   I      T  I  i. LULUS   POISSON'S  RATIO  (VARIABLE)   ) 

IF  (S-SS)  23,  23,  . 
2  2  Z=.b-(  (  .!j-F)  *SS*US)/(  S«C  ) 

00  10  21 
23  Z  =  P 
21  WRITE  16,12)  S,  D,  DS,  DT,  P,  Z 

12  FORMAT  (  ?X,  E10.4,  AX,  3(6X,  E10.4),  2E15.4) 
*RITE  (6,60) 
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60  FORMAT  [10X.14H  COR        ,rY,8Xf22h  BUCKLi      EFFICIENTS,  luX, 
i43H   RIGIDITY  PARAMI  I  CORRECTED  STRESS  /  13X,  96H  PARAMETE 

i,  [BLE  1  IBLE  ELASTIC  COMPRESSIBLE  INCOMPRESSIBLE 

1  CONST.  P   VARIABLE  P   ) 

G=D/(2.*(1 .+P) ) 

GS  =  DS/(  :'••-(  L.-2.  *P)*[  S/l  ) 

GT=DT/ (3. -(1.-2. *P) *DT/I  I 

XK=u/{3.*( L.-2.*P) ) 
-3.*XK 

XH=(4.*   *(   •»,))/(  XN+4.*GS  ) 

YG=  (2.*GS*(  <  \-2. *GS)  )  /  (  ■>'.  *  '•  .  *GS  ) 

S1=.4714*S 
[R=(3.*(GT-GS) )/(S*S*<XN+4.*GS)*(XN+3.*GS+GT) ) 

A11=(2.*GS i •  '  (R*( 2.*GS*S)**2)/XH 

A22=(2.*GS+YG  +  X  -  1  (  IN+i  .*GS  l*S )*»2 )/XH 

A33=4.*bS/XH 

A12=(YG-   -  .*GS*(XN+  . *GS) *S*S )/XH 

B22=l.-.75*( (DS-DT)/DS) 

PH=uS/U 

XE=D/( (1 .-P*P)*> HJ 

A1=.5*(A22*A33+4.*A11*A;        .*A3  -,1  *B*  -4.*A12«A12*B*B+A11*A3 
L3*B**4) 

A2=2.*Bj  HA11*B*B1 
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A3=2.*XE*B**4*(A11+A22-A33-2.*A12) 

A4=.5*B*B*XE*( l.+B*B)*A33 

ao  =  Z  .  *  (  B*  _  -  !E  )  »*2 

B1=.5*(B22+4.*B22*B*B-2.*B*B+B**4)*PH 

B2=2.*B*B*(B22+B*B] 

B3  =  <L . *(B22-J  . ) *B**4 

B4=. 5*B»B*( L.+B*E )*PK 

B5=2  . *B**4 

Cl=( l.+B*B) **2 

C2=( l.+B*B) 

PA  =  2 .»S*{ i.-P*P) 
70    R=0. 

50     IF     (R-1.1J      i     ,4     ,20 
40    XKP1  =     (l.+B*B)*{Al*R+A2)-A3 

XKP2=K*(  Al*  H  A2  )  +A4*R+A5 

XKP=XKPL/XKP; 

YKP1=(1.+B*B)*(B1*R+B2)- 

YKP2=R*PH«{31*R+B2)+B4*  + 

YKP=(PH-      , /YKP2 

YKl=Cl/< R*C;  *B*B  ) 

ZN=.75/( l.-P*P) 

Z1=PA/XKP 

Z3=1.5*S/YI 
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Z2=/N*S 

ZM=. 7b/( L.-Z*Z) 

Z4=ZM*S 

(6,     LJ        .     XKP,     YKP,    YK1,     Zl,  ,  ,     M 

6  1    FORMAT     (15Xf     F6.2,    2X,     3E12.5,     2X,    2E11.3,    4X,  :L.3) 

R=R+.2 

2  0    C0N1  [I 
RETUR  . 

E  . 


B.      MATERIAL  PROPERTIES 

Here,    the  material    properties    are   given   for    17-7   PH  stainless 
steel.     The   stress-strain  curve  was   taken   from  the  book  by  Watter   and 
Lincoln    [31]   and   is   plotted   in  Fig.   6.      From  the  curve,    the  various 
moduli   (i.e.,   Young's,    secant,    and   tangent)  were  measured.      Young's 
modulus   was    found   to  be     28.6   x   10     psi,    and   the  secant   and   tangent 
moduli  are   plotted   in  Fig.    7.     Also,    the   0.02%  yield  stress  was   found 
to  be   68,000  psi.      In  addition,    from  reference    [32],   Poisson's   ratio 
was    taken   to  be    0.28. 
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Fig.   6.      Stress-strain  curve   for   17-7   PH  stainless   steel 
in   compression 
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Fig.    7.      Tangent   and   secant   moduli  for    17-7   PH   stainless    steel 
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took  a  leave  of  absence  from  Westinghouse  to  enroll  in  the  Graduate 
School  of  the  University  of  Florida.  The  University  of  Florida  awarded 
him  with  a  NASA  traineeship  and  graduated  him  with  the  degree  of 
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Master  of  Science  in  Engineering  in  December,  1964;  since  then,  Fred 
has  been  studying  in  the  Department  of  Engineering  Mechanics  at  the 
University  of  Florida  in  the  pursuit  of  a  doctorate  degree. 
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